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These lecture notes provide a concise, rapid and pedagogical introduction to several advanced top- 
ics in contemporary cosmology. The discussion of thermal history of the universe, linear perturbation 
theory, theory of CMBR temperature anisotropics and the inflationary generation of perturbation 
are presented in a manner accessible to someone who has done a first course in cosmology. The 
discussion of dark energy is more research oriented and reflects the personal bias of the author. Con- 
tents: (I) The cosmological paradigm and Friedmann model; (II) Thermal history of the universe; 
(III) Structure formation and linear perturbation theories; (IV) Perturbations in dark matter and 
radiation; (V) Transfer function for matter perturbations; (VI) Temperature anisotropies of CMBR; 
(VII) Generation of initial perturbations from inflation; (VIII) The dark energy. 



Contents 



I. The Cosmological Paradigm and Friedmann mo del 



II. The rmal History of the Universe] 



B^^rimardial^Iu^kosvn 



C. Decoupling of matter and radiation! 



III. Structure Formation and Linear Perturbation Theory! 



IV. Pe rturbations in Dark Matter and Radiation! 

A™ Evolution for A ~5> du\ 



B. Evolution for A <C du in the radiation dominated phase I 



ii " , I 

^^F^^hatio^^r^^e^ma^^r^o^mater^rjhascJ 



D. An alternative description of matter-radiation system! 



V. Transfer Function for matter perturbations! 
VI. Temperature anisotropies of CMBRl 



A. CMBR Temperature Anisotropy: More detailed thcorvj 



P^^csc^jtioiiof^Ao'^m^Da^rori^luidf 



VII. Generation of initial perturbations from inflation I 



VIII. The Dark Energy! 



A. Gravitational Holography! 



C^^Ge^mfitri^al^ualitvji^o^ 



D. Gravity as detector of the vacuum energy! 



Acknowledgement! 
References! 



4 
4 
5 
9 

10 

14 
15 
15 
16 
17 

19 

20 
22 
24 

26 

31 
34 
34 
35 
3G 

37 

37 



I. THE COSMOLOGICAL PARADIGM AND FRIEDMANN MODEL 

Observations show that the universe is fairly homogeneous and isotropic at scales larger than about 150/i _1 Mpc 
where 1 Mpc ~ 3x 10 24 cm and h ss 0.7 is a parameter related to the expansion rate of the universe. The conventional 
— and highly successful — approach to cosmology separates the study of large scale (I > lSOh^ 1 Mpc) dynamics of 
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the universe from the issue of structure formation at smaller scales. The former is modeled by a homogeneous and 
isotropic distribution of energy density; the latter issue is addressed in terms of gravitational instability which will 
amplify the small perturbations in the energy density, leading to the formation of structures like galaxies. In such an 
approach, the expansion of the background universe is described by the metric (We shall use units with with c = 1 
throughout, unless otherwise specified): 

ds 2 = dt 2 - a 2 dx 2 = dt 2 - a 2 (t) [d X 2 + S 2 k (x) (dO 2 + sin 2 9d(/) 2 )] (1) 

with Sk(x) — ( sm X; Xi smn x) f° r the three values of the label k = (1,0,-1). The function a(t) is governed by the 
equations: 



8nGp 



d(pa 3 ) = -pda 3 (2) 



The first one relates expansion rate of the universe to the energy density p and fc = 0,±1 is a parameter which 
characterizes the spatial curvature of the universe. The second equation, when coupled with the equation of state 
V = p(p) which relates the pressure p to the energy density, determines the evolution of energy density p = p(a) in 
terms of the expansion factor of the universe. In particular if p = wp with (at least, approximately) constant w then, 
p cx a~ 3 ( 1+tu ) and (if we further assume k = 0, which is strongly favoured by observations) the first equation in Eq.J2J 
gives a cx 

i 2/[3(l+t«)] < We wiU alg0 often 

use the redshift z(t), defined as (1 + z) = a /a(t) where the subscript zero 
denotes quantities evaluated at the present moment, in a k — universe, we can set oq = 1 by rescaling the spatial 
coordinates. 

It is convenient to measure the energy densities of different components in terms of a critical energy density (p c ) 
required to make k = at the present epoch. (Of course, since k is a constant, it will remain zero at all epochs if it 
is zero at any given moment of time.) From Eq.(J2J), it is clear that p c = 3H 2 /8irG where Ho = (a/a)o — called the 
Hubble constant — is the rate of expansion of the universe at present. Numerically 

3H 2 

Pc = 7T-% = L88/l2 x 10_29 § m cm ~ 3 = 2 - 8 x lO n /i 2 M Mpc" 3 

= 1.1 x 10 4 h 2 eV cm~ 3 = 1.1 x 10~ 5 ft, 2 protons cm" 3 (3) 

The variables fl; = pij p c will give the fractional contribution of different components of the universe (i denoting 
baryons, dark matter, radiation, etc.) to the critical density. Observations then lead to the following results: 

(1) Our universe has 0.98 < £ltot ^ 1.08. The value of O tot can be determined from the angular anisotropy spectrum 
of the cosmic microwave background radiation (CMBR; see Section IVl|) and these observations (combined with the 
reasonable assumption that h > 0.5) show[l| that we live in a universe with critical density, so that k = 0. 

(2) Observations of primordial deuterium produced in big bang nucleosynthesis (which took place when the universe 
was about few minutes in age) as well as the CMBR observations showQ that the total amount of baryons in the 
universe contributes about Qb — (0.024±0.0012)/i -2 . Given the independent observations Q which fix h = 0.72±0.07, 
we conclude that tts — 0.04 — 0.06. These observations take into account all baryons which exist in the universe 
today irrespective of whether they are luminous or not. Combined with previous item we conclude that most of the 
universe is non-baryonic. 

(3) Host of observations related to large scale structure and dynamics (rotation curves of galaxies, estimate of 
cluster masses, gravitational lensing, galaxy surveys ..) all suggestQ that the universe is populated by a non-luminous 
component of matter (dark matter; DM hereafter) made of weakly interacting massive particles which does cluster at 
galactic scales. This component contributes about £Idm — 0.20 — 0.35 and has the simple equation of state pdm ~ 0. 
The second equation in Eq.(J2J), then gives pdm oc a -3 as the universe expands which arises from the evolution of 
number density of particles: p = nine 2 cx n cx a~ 3 . 

(4) Combining the last observation with the first we conclude that there must be (at least) one more component to 
the energy density of the universe contributing about 70% of critical density. Early analysis of several observations 
indicated that this component is unclustered and has negative pressure. This is confirmed dramatically by the 
supernova observations (see Ref. 0]; for a critical look at the current data, see Ref. Q)- The observations suggest 
that the missing component has w = p/p < —0.78 and contributes FIde — 0.60 — 0.75. The simplest choice for 
such dark energy with negative pressure is the cosmological constant which is a term that can be added to Einstein's 
equations. This term acts like a fluid with an equation of state pde = ~Pde\ the second equation in Eq.(|2J), then 
gives pde = constant as universe expands. 

(5) The universe also contains radiation contributing an energy density f2#/i 2 = 2.56 x 10~ 5 today most of which 
is due to photons in the CMBR. The equation of state is pn — (l/3)pi?,; the second equation in Eq.J5J), then gives 
pa oc a~ A . Combining it with the result pji cx T 4 for thermal radiation, it follows that T oc a -1 . Radiation is 
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dynamically irrelevant today but since (pr/ Pdm) oc a -1 it would have been the dominant component when the 
universe was smaller by a factor larger than SIdmI^r — 4 x 10 i fl]jM h 2 . 

(6) Taking all the above observations together, we conclude that our universe has (approximately) Qde — 
0.7, fi^M — 0.26, fis ~ 0.04, f^T? ~ 5 x 10~ 5 . All known observations are consistent with such an — admittedly 
weird — composition for the universe. 

Using pnr oc a~ 3 , pa oc a -4 and pnB^constant we can write Eq.J2J in a convenient dimensionless form as 

2 

+ V(q) = E (4) 



where r = Hot, a = aoq(r), J7nr = 

V(q) 









K 


+ I^dm and 


1 




~2 





dqY 
dr J 



XR 



E=-(l- fit*) • 



(5) 



This equation has the structure of the first integral for motion of a particle with energy E in a potential V(q). For 
models with f2 = O^r + ^-de = 1, we can take E — so that (dq/dr) = y/V(q). Based on the observed composition 
of the universe, we can identify three distinct phases in the evolution of the universe when the temperature is less 
than about 100 GeV. At high redshifts (small q) the universe is radiation dominated and q is independent of the 
other cosmological parameters. Then Eq.Q can be easily integrated to give a(t) oc t 1 / 2 and the temperature of the 



universe decreases as T oc t 1 I 2 . As the universe expands, a time will come when (t = i eq , a a. 



, .q and Z — ^cq: 



say) the matter energy density will be comparable to radiation energy density. For the parameters described above, 
(1 + z eq ) = VInr/Qr ~ 4 x 10 4 n DM h 2 . At lower redshifts, matter will dominate over radiation and we will have 
a oc i 2 / 3 until fairly late when the dark energy density will dominate over non relativistic matter. This occurs at a 
redshift of zde where (1 + zde) = (^de/^nr) 1 ^ 3 - For f^E ~ 0.7, JInr ~ 0.3, this occurs at zde ~ 0.33. In this 
phase, the velocity q changes from being a decreasing function to an increasing function leading to an accelerating 
universe. In addition to these, we believe that the universe probably went through a rapidly expanding, inflationary, 
phase very early when T « 10 14 GeV; we will say more about this in Section lYlII (For a textbook description of these 
and related issues, see e.g. Ref. @.) 

Before we conclude this section, wc will briefly mention some key aspects of the background cosmology described 
by a Friedmann model. 

(a) The metric in Eq.([Q) can be rewritten using the expansion parameter a or the redshift z — (ao/a) -1 — 1 as the 
time coordinate in the form 



ds z 



■(.) if 



a dx 



1 



1 



[H- 2 (z)dz 2 - dx 2 ] 



(6) 



This form clearly shows that the only dynamical content of the metric is encoded in the function H(a) = (a /a). An 
immediate consequence is that any observation which is capable of determining the geometry of the universe can only 
provide — at best — information about this function. 

(b) Since cosmological observations usually use radiation received from distant sources, it is worth reviewing briefly 
the propagation of radiation in the universe. The radial light rays follow a trajectory given by 

r cm {z) = S k (a); a= — [ H~ l {z)dz (7) 



a = — H-\z)dz 
ao Jo 

if the photon is emitted at r cm at the redshift z and received here today. Two other quantities closely related to 
r cm (z) are the luminosity distance, d^, and the angular diameter distance dA- If we receive a flux F from a source 



of luminosity L, then the luminosity distance is defined via the relation F = L/ATrd L 
length I subtends a small angle 8, the angular diameter distance is defined via (/ = 
that: 



). If an object of transverse 
). Simple calculation shows 



d L (z) = a r cm (z)(l + z) = a (l + z)S k {a); d A {z) = a r cm (z)(l + z)~ 



(8) 



(c) As an example of determining the spacetime geometry of the universe from observations, let us consider how 
one can determine a(t) from the observations of the luminosity distance. It is clear from the first equation in Eq. © 
that 



H-\z) = 



1 - 



kd\{z) 
a 2 (l + z) 2 



■ - 1/2 d 


\d L (zY 


d 


\d L (z)] 


dz 


.1 + *. 


dz 


.1 + *. 



(9) 



where the last form is valid for a k — universe. If we determine the form of dL{z) from observations — which can 
be done if we can measure the flux F from a class of sources with known value for luminosity L — then we can use 
this relation to determine the evolutionary history of the universe and thus the dynamics. 
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II. THERMAL HISTORY OF THE UNIVERSE 

Let us next consider some key events in the evolutionary history of our universe Q . The most well understood phase 
of the universe occurs when the temperature is less than about 10 12 K. Above this temperature, thermal production of 
baryons and their strong interaction is significant and somewhat difficult to model. We can ignore such complications 
at lower temperatures and — as we shall see — several interesting physical phenomena did take place during the later 
epochs with T < 10 12 . 

The first thing we need to do is to determine the composition of the universe when T w 10 12 K. We will certainly 
have, at this time, copious amount of photons and all species of neutrinos and antineutrinos. In addition, neutrons and 
protons must exist at this time since there is no way they could be produced later on. (This implies that phenomena 
which took place at higher temperatures should have left a small excess of baryons over anti baryons; we do not 
quite understand how this happened and will just take it as an initial condition.) Since the rest mass of electrons 
correspond to a much lower temperature (about 0.5 x 10 10 K), there will be large number of electrons and positrons at 
this temperature but in order to maintain charge neutrality, we need to have a slight excess of electrons over positrons 
(by about 1 part in 10 9 ) with the net negative charge compensating the positive charge contributed by protons. 

An elementary calculation using the known interaction rates show that all these particles are in thermal equilibrium 
at this epoch. Hence standard rules of statistical mechanics allows us to determine the number density (n), energy 
density (p) and the pressure (p) in terms of the distribution function /: 

, , o g f°° (E 2 -m 2 ) 1/2 EdE , . 



Ef(k)d i k=^l ^ ' " I" (11) 

2tH J m exp[(E - fi)/T] ± 1 



(E 2 - m 2 ) 1 / 2 E 2 dE 



P = -J dV(k)Mk) - J 3^/W*k = ^ I \ liE ^ m±l (12) 

Next, we can argue that the chemical potentials for electrons, positrons and neutrinos can be taken to be zero. For 
example, conservation of chemical potential in the reaction e + e~ — > 2j implies that the chemical potentials of electrons 
and positrons must differ in a sign. But since the number densities of electrons and positrons, which arc determined 
by the chemical potential, are very close to each other, the chemical potentials of electrons and positrons must be 
(very closely) equal to each other. Hence both must be (very close to) zero. Similar reasoning based on lepton number 
shows that neutrinos should also have zero chemical potential. Given this, one can evaluate the integrals for all the 
relativistic species and we obtain for the total energy density 

ptotai = e 9i (io) T * + \ g% (io) T i = 5totai (io) 714 ^ 



boson i— fer 



mion 



where 



£total 



E 9b + J2 h F - (14) 



boson fcrmion 



The corresponding entropy density is given by 

— (p + p) = — 
T KH FJ 45 



s - ^(p + p) = -jrvT 3 '' 1 = 9*otai = E 9B + 5 E 9F - ( 15 ) 



boson fcrmion 



A. Neutrino background 

As a simple application of the above result, let us consider the fate of neutrinos in the expanding universe. From the 
standard weak interaction theory, one can compute the reaction rate T of the neutrinos with the rest of the species. 
When this reaction rate fall below the expansion rate H of the universe, the reactions cannot keep the neutrinos 
coupled to the rest of the matter. A simple calculation || shows that the relevant ratio is given by 

r / t \ 3 _ / t 

H~\TAMeVj ~ VL6x 10 10 K 
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Thus, for T < 1.6 x 10 10 K, the neutrinos decouple from matter. At slightly lower temperature, the electrons and 
positrons annihilate increasing the number density of photons. Neutrinos do not get any share of this energy since 
they have already decoupled from the rest of the matter. As a result, the photon temperature goes up with respect 
to the neutrino temperature once the e + e~ annihilation is complete. This increase in the temperature is easy to 
calculate. As far as the photons are concerned, the increase in the temperature is essentially due to the change in the 
degrees of freedom g and is given by: 

(^ 7 ) after = ^before = 1(2 + 2) + 2 = 11 
(^before fafter 2 4 ' 

(In the numerator, one 2 is for electron; one 2 is for positron; the 7/8 factor arises because these are fermions. The 
final 2 is for photons. In the denominator, there are only photons to take care of.) Therefore 

n xl/3 x nX l/3 



(aT 7 ) aftcI . — ( —j (aT 7 )before = y~ ) ( a ^w)before 
U XV3 



—J (ar i/ ) a ft er ~1.4(aT y ) after . (18) 

The first equality is from Eq. Q17|): the second arises because the photons and neutrinos had the same temperature 
originally; the third equality is from the fact that for decoupled neutrinos aT u is a constant. This result leads to the 
prediction that, at present, the universe will contain a bath of neutrinos which has temperature that is (predictably) 
lower than that of CMBR. The future detection of such a cosmic neutrino background will allow us to probe the 
universe at its earliest epochs. 



B. Primordial Nucleosynthesis 

When the temperature of the universe is higher than the binding energy of the nuclei (~ MeV) , none of the heavy 
elements (helium and the metals) could have existed in the universe. The binding energies of the first four light nuclei, 
2 H, 3 H, 3 He and 4 He are 2.22 MeV, 6.92 MeV, 7.72 MeV and 28.3 MeV respectively. This would suggest that these 
nuclei could be formed when the temperature of the universe is in the range of (1 — 30)MeV. The actual synthesis 
takes place only at a much lower temperature, T nuc = T n ~ O.lMeV. The main reason for this delay is the 'high 
entropy' of our universe, i.e., the high value for the photon-to-baryon ratio, ?7 _1 . Numerically, 

' = 5 - 5 * 10 "° (w) ; n »*' - M5 * 10 " {^k) 3 '"« < w > 

To see this, let us assume, for a moment, that the nuclear (and other) reactions are fast enough to maintain thermal 
equilibrium between various species of particles and nuclei. In thermal equilibrium, the number density of a nuclear 
species A Nz with atomic mass A and charge Z will be 



'm A T\ 
n A = 9 A 1 —~ — I exp 



V T 



From this one can obtain the equation for the temperature T A at which the mass fraction of a particular species-A 
will be of order unity (Xa — 1). We find that 

B A /(A~1) 

A ~ \n(r ] - 1 ) + 1.5 hx(m B /T) { ' 

where Ba is the binding energy of the species. This temperature will be fairly lower than Ba because of the large 
value of 77 _1 . For 2 H, 3 He and 4 He the value of Ta is 0.07MeV, O.llMeV and 0.28MeV respectively. Comparison 
with the binding energy of these nuclei shows that these values are lower than the corresponding binding energies Ba 
by a factor of about 10, at least. 

Thus, even when the thermal equilibrium is maintained, significant synthesis of nuclei can occur only at T < 0.3McV 
and not at higher temperatures. If such is the case, then we would expect significant production (Xa < 1) of nuclear 
species-A at temperatures T < Ta- It turns out, however, that the rate of nuclear reactions is not high enough to 
maintain thermal equilibrium between various species. We have to determine the temperatures up to which thermal 
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equilibrium can be maintained and redo the calculations to find non-equilibrium mass fractions. The general procedure 
for studying non equilibrium abundances in an expanding universe is based on rate equations. Since we will require 
this formalism again in Section lll CI ffor the study of recombination), we will develop it in a somewhat general context. 

Consider a reaction in which two particles 1 and 2 interact to form two other particles 3 and 4. For example, 
n+ v e ^ p + e constitutes one such reaction which converts neutrons into protons in the forward direction and protons 
into neutrons in the reverse direction; another example we will come across in the next section isp + e^H +7 
where the forward reaction describes recombination of electron and proton forming a neutral hydrogen atom (with 
the emission of a photon), while the reverse reaction is the photoionisation of a hydrogen atom. In general, we are 
interested in how the number density n\ of particle species 1, say, changes due to a reaction of the form 1 + 2^3 + 4. 

We first note that even if there is no reaction, the number density will change as n\ oc a~ 3 due to the expansion 
of the universe; so what we are really after is the change in n^o? . Further, the forward reaction will be proportional 
to the product of the number densities n\ni while the reverse reaction will be proportional to 713114. Hence we can 
write an equation for the rate of change of particle species n\ as 

1 d(n ia 3 ) 

—- pyAn^n^ — n\n2). (22) 



a 3 dt 

The left hand side is the relevant rate of change over and above that due to the expansion of the universe; on the 
right hand side, the two proportionality constants have been written as fi and (Ap), both of which, of course, will 
be functions of time. (The quantity p has the dimensions of cm 3 s _1 , so that np has the dimensions of s ; usually 
fi ~ crv where a is the cross-section for the relevant process and v is the relative velocity.) The left hand side has to 
vanish when the system is in thermal equilibrium with m = n° q , where the superscript 'eq' denotes the equilibrium 
densities for the different species labeled by i = 1 — 4. This condition allows us to rewrite A as A = n^n^ 1 / (n^n^). 
Hence the rate equation becomes 

1 d(nia 3 ) _ cq cq / n^tii nin 2 

Q 77 t m l n 2 I oq cq cq cq 

a 6 dt \na H n A n-, n. 



fiU 1 U 2 I oq oq oq cq ) • (23) 



In the left hand side, one can write [d/dt] = Ha(d/da) which shows that the relevant time scale governing the 
process is H~ 1 . Clearly, when H/nfi ^> 1 the right hand side becomes ineffective because of the (p/H) factor and 
the number of particles of species 1 does not change. We see that when the expansion rate of the universe is large 
compared to the reaction rate, the given reaction is ineffective in changing the number of particles. This certainly does 
not mean that the reactions have reached thermal equilibrium and rn — n- q ; in fact, it means exactly the opposite: 
The reactions are not fast enough to drive the number densities towards equilibrium densities and the number densities 
"freeze out" at non equilibrium values. Of course, the right hand side will also vanish when ni = n° q which is the 
other extreme limit of thermal equilibrium. 

Having taken care of the general formalism, let us now apply it to the process of nucleosynthesis which requires 
protons and neutrons combining together to form bound nuclei of heavier elements like deuterium, helium etc.. The 
abundance of these elements are going to be determined by the relative abundance of neutrons and protons in the 
universe. Therefore, we need to first worry about the maintenance of thermal equilibrium between protons and the 
neutrons in the early universe. As long as the inter-conversion between n and p through the weak interaction processes 
(y + n <-* p + e), (e + n <-> p + 77) and the 'decay' (n <-> p + e + 77), is rapid (compared to the expansion rate of the 
universe), thermal equilibrium will be maintained. Then the equilibrium (n/p) ratio will be 

= ^=exp(-Q/T), (24) 

Aj, 

where Q = m„ — m p = 1.293 MeV. At high (T 3> Q) temperatures, there will be equal number of neutrons and 
protons but as the temperature drops below about 1.3 MeV, the neutron fraction will start dropping exponentially 
provided thermal equilibrium is still maintained. To check whether thermal equilibrium is indeed maintained, we 
need to compare the expansion rate with the reaction rate. The expansion rate is given by H = (SnGp/S) 1 ^ 2 where 
p = (it 2 /30)gT 4 with g i=s 10.75 representing the effective relativistic degrees of freedom present at these temperatures. 
At T = Q, this gives ff» 1.1 s , The reaction rate needs to be computed from weak interaction theory. The neutron 
to proton conversion rate, for example, is well approximated by 



A„ p w 0.29 s _1 



Q 



T \T 



(25) 



At T = Q, this gives A w 5 s , slightly more rapid than the expansion rate. But as T drops below Q, this decreases 
rapidly and the reaction ceases to be fast enough to maintain thermal equilibrium. Hence wc need to work out the 
neutron abundance by using Eq. (|23|l . 
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FIG. 1: The evolution of mass fraction of different species during nucleosynthesis 



Using ti\ — n n , n 3 = n p and n 2 , n 4 = n; where the subscript / stands for the leptons, Eq. I|23[l becomes 

1 d{n n a 3 ) ( n p n^ \ 

"= M«j I - JL w~ ~ n n ■ (26) 



a 3 dt 

We now use Eq. (|24H . write (n^ q /i) = X np which is the rate for neutron to proton conversion and introduce the 
fractional abundance X n = n n /(n n + rip). Simple manipulation then leads to the equation 



dX n 
~dt 



X np ((1 - X n )e~ Q / T - V„) . (27) 



Converting from the variable t to the variable s — (Q/T) and using (d/dt) = —HT(d/dT), the equations we need to 
solve reduce to 

-Hs^=\ np ((l-X n )e- S -X n ); H = (1.1 sec" 1 ) s" 4 ; X np = °' 29 J - [s 2 + 6s + 12] . (28) 

It is now straightforward to integrate these equations numerically and determine how the neutron abundance changes 
with time. The neutron fraction falls out of equilibrium when temperatures drop below 1 MeV and it freezes to about 
0.15 at temperatures below 0.5 MeV. 

As the temperature decreases further, the neutron decay with a half life of r n ~ 886.7 sec (which is not included in 
the above analysis) becomes important and starts depicting the neutron number density. The only way neutrons can 
survive is through the synthesis of light elements. As the temperature falls further to T — Thc — 0.28MeV, significant 
amount of He could have been produced if the nuclear reaction rates were high enough. The possible reactions which 
produces 4 He are [D(D,n) 3 He{D,p) 4 He, D{D,p) 3 H{D,n) 4 He,D{D,j) 4 He\. These are all based on D, 3 He and 
3 H and do not occur rapidly enough because the mass fraction of D, 3 He and 3 H are still quite small [10~ 12 , 10 -19 
and 5 x 10~ 19 respectively] at T ~ 0.3MeV. The reactions n + p ^ d + -f will lead to an equilibrium abundance ratio 
of deuterium given by 



n p n n _ 4 / m p m n \ 3/2 (27rfc B T) 3/2 ^ B/kBT _ ^ 
n d n 3 \ md J (2irhf n 



25.82 - In n B h 2 T^ 2 



2.58 
IT 



in 



(29) 



The equilibrium deuterium abundance passes through unity (for fish 2 = 0.02) at the temperature of about 0.07 MeV 
which is when the nucleosynthesis can really begin. 
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So we need to determine the neutron fraction at T = 0.07 MeV given that it was about 0.15 at 0.5 MeV. During 
this epoch, the time-temperature relationship is given by t — 130 sec (T/0.1 McV)~ 2 . The neutron decay factor 
is exp(— t/r n ) f» 0.74 for T = 0.07 MeV. This decreases the neutron fraction to 0.15 x 0.74 = 0.11 at the time of 
nucleosynthesis. When the temperature becomes T < 0.07MeV, the abundance of D and 3 H builds up and these 
elements further react to form 4 He. A good fraction of D and 3 H is converted to 4 He (See Fig^ which shows the 
growth of deuterium and its subsequent fall when helium is built up). The resultant abundance of 4 He can be easily 
calculated by assuming that almost all neutrons end up in 4 He. Since each 4 He nucleus has two neutrons, (n n /2) 
helium nuclei can be formed (per unit volume) if the number density of neutrons is n n . Thus the mass fraction of 
4 He will be 

F=^Zg)= 2{n / P) =2 Xc (30) 
n n + n p l + (n/p) 

where x c — n/(n + p) is the neutron abundance at the time of production of deuterium. For VLsh 2 = 0.02, x c « 0.11 
giving Y » 0.22. Increasing baryon density to fish 2 = 1 will make Y » 0.25. An accurate fitting formula for the 
dependence of helium abundance on various parameters is given by 

Y = 0.226 + 0.025 log 7710 + 0.0075( ff * - 10.75) + 0.014(r 1/2 (n) - 10.3 min) (31) 

where 7710 measures the baryon-photon ratio today via Eq. (|19(l and is the effective number of relativistic degrees of 
freedom contributing to the energy density and 1^/2 (77) is the neutron half life. The results (of a more exact treatment) 
are shown in Fig. ^ 

As the reactions converting D and 3 H to 4 He proceed, the number density of D and 3 H is depleted and the reaction 
rates - which are proportional to T oc Xa^u^) < av > - become small. These reactions soon freeze-out leaving a 
residual fraction of D and 3 H (a fraction of about 10 -5 to 10~ 4 ). Since r oc rj it is clear that the fraction of (D 3 H) 
left unreacted will decrease with 77. In contrast, the 4 He synthesis - which is not limited by any reaction rate - is 
fairly independent of 77 and depends only on the (n/p) ratio at T ~ O.lMeV. The best fits, with typical errors, to 
deuterium abundance calculated from the theory, for the range 77 = (10~ 10 — 10~ 9 ) is given by 



>L. S ||) =3.6X10— -L-JLJ-" (32, 



The production of still heavier elements - even those like 16 C, le O which have higher binding energies than 4 He 
- is suppressed in the early universe. Two factors are responsible for this suppression: (1) For nuclear reactions to 
proceed, the participating nuclei must overcome their Coulomb repulsion. The probability to tunnel through the 
Coulomb barrier is governed by the factor F = cxp[— 2A 1 / 3 (ZiZ 2 ) 2/3 (T/lMcV)- 1 / 3 } where A^ 1 = A\ x + A^ 1 . For 
heavier nuclei (with larger Z), this factor suppresses the reaction rate. (2) Reaction between helium and proton would 
have led to an element with atomic mass 5 while the reaction of two helium nuclei would have led to an element with 
atomic mass 8. However, there are no stable elements in the periodic table with the atomic mass of 5 or 8! The 8 Be, 
for example, has a half life of only 10 -16 seconds. One can combine 4 He with 8 Be to produce 12 C but this can occur at 
significant rate only if it is a resonance reaction. That is, there should exist an excited state 12 C nuclei which has an 
energy close to the interaction energy of 4 He + 8 Be. Stars, incidentally, use this route to synthesize heavier elements. 
It is this triple-alpha reaction which allows the synthesis of heavier elements in stars but it is not fast enough in the 
early universe. (You must thank your stars that there is no such resonance in 16 or in 20 Ne — which is equally 
important for the survival of carbon and oxygen.) 

The current observations indicate, with reasonable certainty that: (i) (D/H) > 1 x 10~ 5 . (ii) [(D + 3 He)/H] 
~ (1 — 8) x 10~ 5 and (iii) 0.236 < ( 4 He/H) < 0.254. These observations are consistent with the predictions if 
10.3 min < r < 10.7 min, and rj = (3 - 10) x 10~ 10 . Using rj = 2.68 x 10 -8 fJ B /i 2 , this leads to the important 
conclusion: 0.011 < 5I73/7 2 < 0.037. When combined with the broad bounds on h, 0.6 < h < 0.8, say, we can constrain 
the baryonic density of the universe to be: 0.01 < < 0.06. These are the typical bounds on fig available today. It 
shows that, if fitotai — 1 then most of the matter in the universe must be non baryonic. 

Since the 4 He production depends on g, the observed value of 4 He restricts the total energy density present at 
the time of nucleosynthesis. In particular, it constrains the number (N u ) of light neutrinos (that is, neutrinos with 
i; lMeV which would have been relativistic at T ~ lMeV). The observed abundance is best explained by N u = 3, 
is barely consistent with N u = 4 and rules out N v > 4. The laboratory bound on the total number of particles 
including neutrinos, which couples to the Z° boson is determined by measuring the decay width of the particle Z°; 
each particle with mass less than (m z /2) ~ 46 GeV contributes about 180 MeV to this decay width. This bound is 
N v = 2.79 ± 0.63 which is consistent with the cosmological observations. 
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C. Decoupling of matter and radiation 



In the early hot phase, the radiation will be in thermal equilibrium with matter; as the universe cools below 
kgT ~ (e Q /10) where e a is the binding energy of atoms, the electrons and ions will combine to form neutral atoms 
and radiation will decouple from matter. This occurs at T^ cc ~ 3 x 10 3 K. As the universe expands further, these 
photons will continue to exist without any further interaction. It will retain thermal spectrum since the redshift of the 
frequency v cx a -1 is equivalent to changing the temperature in the spectrum by the scaling T cx (1/a). It turns out 
that the major component of the extra-galactic background light (EBL) which exists today is in the microwave band 
and can be fitted very accurately by a thermal spectrum at a temperature of about 2.73 K. It seems reasonable to 
interpret this radiation as a relic arising from the early, hot, phase of the evolving universe. This relic radiation, called 
cosmic microwave background radiation, turns out to be a gold mine of cosmological information and is extensively 
investigated in recent times. We shall now discuss some details related to the formation of neutral atoms and the 
decoupling of photons. 

The relevant reaction is, of course, e + p ^ H + 7 and if the rate of this reaction is faster than the expansion rate, 
then one can calculate the neutral fraction using Saha's equation. Introducing the fractional ionisation, X,-, for each 
of the particle species and using the facts n p = n e and n p + n B = rig, it follows that X p = X e and Xh = {tih/tib) 
= 1 — X e . Saha's equation now gives 

i— 5-^ £* 3.8477(T/m e ) 3/2 exp(B/T) (33) 

where r\ = 2.68 x 10 _8 (ris/i 2 ) is the baryon-to-photon ratio. We may define T a tom as the temperature at which 90 
percent of the electrons, say, have combined with protons: i.e. when X e = 0.1. This leads to the condition: 

{VtBh 2 )- 1 ^ exp [-13.6T- 1 ] = 3.13 x 10" 18 (34) 

where r = (T/leV). For a given value of (fish 2 ), this equation can be easily solved by iteration. Taking logarithms 
and iterating once we find r _1 = 3.084 — 0.0735 ln(ils/i 2 ) with the corresponding redshift (1 + z) = (T/Tq) given by 

(1 + z) = 1367[1 - 0.024 hx{fl B h 2 )]- x . (35) 

For Q B h 2 = 1,0.1,0.01 we get T atom = 0.324eV, 0.307eV, 0.292eV respectively. These values correspond to the 
redshifts of 1367, 1296 and 1232. 

Because the preceding analysis was based on equilibrium densities, it is important to check that the rate of the 
reactions p + e<-^H + r yis fast enough to maintain equilibrium. For fish 2 « 0.02, the equilibrium condition is only 
marginally satisfied, making this analysis suspect. More importantly, the direct recombination to the ground state of 
the hydrogen atom — which was used in deriving the Saha's equation — is not very effective in producing neutral 
hydrogen in the early universe. The problem is that each such recombination releases a photon of energy 13.6 eV 
which will end up ionizing another neutral hydrogen atom which has been formed earlier. As a result, the direct 
recombination to the ground state does not change the neutral hydrogen fraction at the lowest order. Recombination 
through the excited states of hydrogen is more effective since such a recombination ends up emitting more than one 
photon each of which has an energy less than 13.6 eV. Given these facts, it is necessary to once again use the rate 
equation developed in the previous section to track the evolution of ionisation fraction. 

A simple procedure for doing this, which captures the essential physics, is as follows: We again begin with Eq. ill'MI) 
and repeating the analysis done in the last section, now with n\ = n e ,n2 = n p ,U3 — riR- and m — n 7 , and defining 
X e — n e / (n e + uh) = n p /nn one can easily derive the rate equation for this case: 

^ = [(3(1 - X e ) - an b X 2 ] = a ^(1 - X e ) - n b X 2 ^j . (36) 

This equation is analogous to Eq. (|27|l : the first term gives the photoionisation rate which produces the free electrons 
and the second term is the recombination rate which converts free electrons into hydrogen atom and we have used the 
fact n e = n b X e etc.. Since we know that direct recombination to the ground state is not effective, the recombination 
rate a is the rate for capture of electron by a proton forming an excited state of hydrogen. To a good approximation, 
this rate is given by 

a = 9.78r 2 c(pj 1/ \n(?) (37) 
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where ro = e 2 /m e c 2 is the classical electron radius. To integrate Eq. I| .'Still we also need to know j3/a. This is easy 
because in thermal equilibrium the right hand side of Eq. I|36[) should vanish and Saha's equation tells us the value of 
X e in thermal equilibrium. On using Eq. (|33[1 . this gives 



a 



2vr 



3/2 



exp[-(B/T)]. 



(38) 



We can now integrate Eq. I|36[l using the variable B/T just as we used the variable Q/T in solving Eq. 1)2 7[). The 
result shows that the actual recombination proceeds more slowly compared to that predicted by the Saha's equation. 
The actual fractional ionisation is higher than the value predicted by Saha's equation at temperatures below about 
1300. For example, at z = 1300, these values differ by a factor 3; at z ~ 900, they differ by a factor of 200. The 
value of T at om, however, does not change significantly. A more rigorous analysis shows that, in the redshift range of 
800 < z < 1200, the fractional ionisation varies rapidly and is given (approximately) by the formula, 



X e = 2.4 x 10" 



; (^NR.fe 2 ) 1/2 



dooo) 



12.75 



(39) 



This is obtained by fitting a curve to the numerical solution. 

The formation of neutral atoms makes the photons decouple from the matter. The redshift for decoupling can be 
determined as the epoch at which the optical depth for photons is unity. Using Eq. (|39|l , we can compute the optical 
depth for photons to be 



n(t)X e (t)a T dt 



nl : I A', ( : )o-i ( ^ 



dz ~ 0.37 



1000 



14.25 



(40) 



where we have used the relation H$dt = — f2 NR z~~ 5 > dz which is valid for z ^> 1. This optical depth is unity at 
Zdcc = 1072. From the optical depth, we can also compute the probability that the photon was last scattered in the 
interval [z, z + dz). This is given by (exp — r) (dr j dz) which can be expressed as 



dr 

P(z) = e~ T ^ = 5.26 x 10~ 3 
dz 



dooo) 



13.25 



exp 



-0.37 



1000 



14.25 



(41) 



This P(z) has a sharp maximum at z ~ 1067 and a width of about Az = 80. It is therefore reasonable to assume 
that decoupling occurred at z ~ 1070 in an interval of about Az ~ 80. We shall see later that the finite thickness of 
the surface of last scattering has important observational consequences. 



III. STRUCTURE FORMATION AND LINEAR PERTURBATION THEORY 



Having discussed the evolution of the background universe, we now turn to the study of structure formation. Before 
discussing the details, let us briefly summarise the broad picture and give references to some of the topics that we will 
not discuss. The key idea is that if there existed small fluctuations in the energy density in the early universe, then 
gravitational instability can amplify them in a well-understood manner leading to structures like galaxies etc. today. 
The most popular model for generating these fluctuations is based on the idea that if the very early universe went 
through an inflationary phase 0], then the quantum fluctuations of the field driving the inflation can lead to energy 
density fluctuations |l0L Hi} . It is possible to construct models of inflation such that these fluctuations are described 
by a Gaussian random field and are characterized by a power spectrum of the form P(k) = Ak n with n ~ 1 (see Sec. 
IVIIl! . The models cannot predict the value of the amplitude A in an unambiguous manner but it can be determined 
from CMBR observations. The CMBR observations are consistent with the inflationary model for the generation of 
perturbations and gives A ~ (28.3/i~ 1 Mpc) 4 and n = 0.97 ± 0.023 (The first results were from COBE [±2 and WMAP 
has reconfirmed them with far greater accuracy). When the perturbation is small, one can use well defined linear 
perturbation theory to study its growth. But when S f=a (Sp/p) is comparable to unity the perturbation theory breaks 
down. Since there is more power at small scales, smaller scales go non-linear first and structure forms hierarchically. 
The non linear evolution of the dark matter halos (which is an example of statistical mechanics of self gravitating 
systems; see e.g.0) can be understood by simulations as well as theoretical models based on approximate ansatz 
|14| and nonlinear scaling relations [l5j . The baryons in the halo will cool and undergo collapse in a fairly complex 
manner because of gas dynamical processes. It seems unlikely that the baryonic collapse and galaxy formation can be 
understood by analytic approximations; one needs to do high resolution computer simulations to make any progress 
[l6j| . All these results are broadly consistent with observations. 
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As long as these fluctuations are small, one can study their evolution by linear perturbation theory, which is 
what we will start with |17|. The basic idea of linear perturbation theory is well defined and simple. We perturb 
the background FRW metric by gf RW — > gf RW + ^ifc an d a ls° perturb the source energy momentum tensor by 



rpFRW 



„FRW 
9ik 

T^, RW + STik. Linearising the Einstein's equations, one can relate the perturbed quantities by a relation 



of the form C(gf k RW )hn : = STik where C is second order linear differential operator depending on the back ground 
metric gf k RW ■ Since the background is maximally symmetric, one can separate out time and space; for e.g, if k = 0, 
simple Fourier modes can be used for this purpose and we can write down the equation for any given mode, labelled 
by a wave vector k as: 



C(a(t),k)h ab {t,k) =8T ab (t,k) 

To every mode we can associate a wavelength normalized to today's value: A(t) 
mass scale which is invariant under expansion: 



(42) 

(27r/fc)(l + z) _1 and a corresponding 



M = 



X(t) 



47rp / A, 



^ ) =1.5 x 10 n M Q (fi m /j 2 ) 



IMpc 



(43) 



The behaviour of the mode depends on the relative value of X(t) as compared to the Hubble radius dn(t) = (a/a) -1 . 
Since the Hubble radius: dn{t) oc t while the wavelength of the mode: X(t) cx a(t) cx (i 1 / 2 ,/; 2 / 3 ) in the radiation 
dominated and matter dominated phases it follows that X(t) > dn(t) at sufficiently early times. When X(t) ~ dn(t), 
we say that the mode is entering the Hubble radius. Since the Hubble radius at z — z eq is 




= 14Mpc(ft NR ft, 2 )- 



(44) 



it follows that modes with Ao > A cq enter Hubble radius in MD phase while the more relevant modes with A < A oq 
enter in the RD phase. Thus, for a given mode we can identify three distinct phases: First, very early on, when 
A > d 



H 



z > z 



eg 



the dynamics is described by general relativity. In this stage, the universe is radiation dominated, 
gravity is the only relevant force and the perturbations are linear. Next, when A < d# and z > z eq one can describe 
the dynamics by Newtonian considerations. The perturbations are still linear and the universe is radiation dominated. 
Finally, when A < d# , z < z eq we have a matter dominated universe in which we can use the Newtonian formalism; 
but at this stage — when most astrophysical structures form — we need to grapple with nonlinear astrophysical 
processes. 

Let us now consider the metric perturbation in greater detail. When the metric is perturbed to the form: g a b — > 
gab + h a b the perturbation can be split as h ab = (hoo, ho a = w ai hap). We also know that any 3- vector w(x) can be 
split as w = + w" in which w" = V$" is curl-free (and carries one degree of freedom) while is divergence- free 
(and has 2 degrees of freedom). This result is obvious in k— space since we can write any vector w(k) as a sum of two 
terms, one along k and one transverse to k: 



w(k) = w^k) + w ± (k) = k 



k 2 



,, , ,i k- w(k) 

W ( k ) - k I 1,2 



k X 



,11 



0; k • w- 1 = 



(45) 



along k 



transverse to k 



Fourier transforming back, we can split w into a curl- free and divergence-free parts. Similar decomposition works for 
h a p by essentially repeating the above analysis on each index. We can write: 



1 



Kp = il>5 a p + (V a uj + Vpui) + - -S a pV 2 $i ^1 + 2 + 1 + 2 = 6 



(46) 



trace tracelcss from vector 



tracclcss from scalar 



The 11^ is divergence free and h^p is traceless and divergence free. Thus the most general perturbation h ab (ten 
degrees of freedom) can be built out of 



hab = ("no, V* = w a ,h a p) = [/loo^'^w- 1 ),^,*!,^,/^)] =>• [1,(1,2), (1,1,2,2)] 



(47) 



We now use the freedom available in the choice of four coordinate transformations to set four conditions: <P = $i = 
and = thereby leaving six degrees of freedom in («oo = 2$,^»,w , fyja") as nonzero. Then the perturbed line 
element takes the form: 

ds 2 = a 2 (r/) [{1 + 2$(x,?7)} d V 2 - 2w^(x,ry) d V dx a - {(1 - 2^(x,r?)) 6 a p + 2^(x,ry)} dx a dx } (48) 
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To make further simplification we need to use two facts from Einstein's equations. It turns out that the Einstein's 
equations for w and h^p decouple from those for (&,ifj). Further, in the absence of anisotropic stress, one of the 
equations give ip = If we use these two facts, we can simplify the structure of perturbed metric drastically. As far 
as the growth of matter perturbations are concerned, we can ignore and hzg and work with a simple metric: 



ds 2 = a 2 (i)[(l + 2$)(V - (1 - 2<S>)8 a pdx a dx 13 ] 



(49) 



with just one perturbed scalar degree of freedom in $. This is what we will study. 

Having decided on the gauge, let us consider the evolution equations for the perturbations. While one can directly 
work with the Einstein's equations, it turns out to be convenient to use the equations of motion for matter variables, 
since we are eventually interested in the matter perturbations. In what follows, we will use the over-dot to denote 
(d/di]) so that the standard Hubble parameter is H = (\/a)(da/dt) = a/a 2 . With this notation, the continuity 
equation becomes: 



p + 3[aH-<f>)(p + p) = -V a [(p+p)v a ] (50) 

Since the momentum flux in the relativistic case is (p + p)v a , all the terms in the above equation are intuitively 
obvious, except probably the $ term. To see the physical origin of this term, note that the perturbation in Eq. I|49|) 
changes the factor in front of the spatial metric from a 2 to et 2 (l — 2$) so that In a — > In a — $; hence the effective 
Hubble parameter from (a/a) to (a/a) — $ which explains the extra $ term. This is, of course, the exact equation for 
matter variables in the perturbed metric given by Eq. Q49[) : but we only need terms which are of linear order. Writing 
the curl-free velocity part as v a — V Q w, the linearised equations, for dark matter (with p = 0) and radiation (with 
p = (1/3)/?) perturbations are given by: 



8m. 



'L 

dr\ 



3$: 



3 • _ d_ ( 6n R ^ _ 

4 drj \n R 



= V 2 v R + 3$ 



(51) 



where n m and ur are the number densities of dark matter particles and radiation. The same equations in Fourier 
space [using the same symbols for, say, 8(t, x) or 8(t, k)] are simpler to handle: 



k d 
ar\ 

Note that these equations imply 



= —k v„ 



3$: 



_ d_ ( 6n R \ _ 
drj\n R J 



drj 



Sn r< 



d_ 

drj 



d 

dq 



In I — 



k 2 v R + 3$ 



-k 2 (v R - v m ) 



(52) 



(53) 



For long wavelength perturbations (in the limit of k —> 0), this will lead to the conservation of perturbation S(s/n m ) 
in the entropy per particle. 

Let us next consider the Euler equation which has the general form: 



d v [( P + p)v a ] = -( P + p)V a <S>-V a p~4aH(p + p)v c 



(54) 



Once again each of the terms is simple to interpret. The (p+p) arises because the pressure also contributes to inertia 
in a relativistic theory and the factor 4 in the last term on the right hand side arises because the term v a d v (p + p) 
on the left hand side needs to be compensated. Taking the linearised limit of this equation, for dark matter and 
radiation, we get: 



= $ — aHv n 



VR = ^+ ~Sr 



(55) 



Thus we now have four equations in Eqs. (|52fl . (|55|l for the five variables (S m , Sr, v m , vr, $). All we need to do is to 
pick one more from Einstein's equations to complete the set. The Einstein's equations for our perturbed metric are: 



o component : fc 2 $ + 3- ( 4> + j = -AirGa 2 ^ p a 5a = -^Ga 2 p bg 5 total 



° component : $ + -$ = -inGa 2 + p) 



aVa; 



\7v 



" component : 3-$ + 2-$ ^$ + $ = 4irGa 2 8p 



(56) 
(57) 
(58) 
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where A denotes different components like dark matter, radiation etc. Using Eq. (|57|l in Eq. (|56(l we can get a modified 
Poisson equation which is purely algebraic: 

-fc 2 $ = AnGa 2 (paSa - 3 Q (fiA + Pa)v a ^ (59) 

which once again emphasizes the fact that in the relativistic theory, both pressure and density act as source of gravity. 

To get a feel for the solutions let us consider a flat universe dominated by a single component of matter with the 
equation of state p = wp. (A purely radiation dominated universe, for example, will have w = 1/3.) In this case the 
Friedmann background equation gives p oc a^ 3 ^ 1+w ^ and 

£ = 2 £ = 2(1 - 3 W ) 

a (l + 3w)??' a (l + 3iu)V 1 ' 

The equation for the potential <& can be reduced to the form: 

^ + 6(l + W )g + fc2w(& = 
1 + 3u> ry 

The second term is the damping due to the expansion while last term is the pressure support that will lead to 
oscillations. Clearly, the factor kr\ determines which of these two terms dominates. When the pressure term dominates 
(krj » 1), we expect oscillatory behaviour while when the background expansion dominates (krj <C 1), we expect the 
growth to be suppressed. This is precisely what happens. The exact solution is given in terms of the Bcsscl functions 



<7i(k)J„ /2 (V^?) + C 2 (k)n /2 (V^M 5 + 3w 

= ' ^72 > V = 1~Xq ( 62 ) 

rfl 2 - 1 + 3w 

From the theory of Besscl functions, we know that: 

L m o^ ( ^ ^r( y /2 + i) ; }^ Y »^ x) *-^ (63) 

This shows that if we want a finite value for $ as r/ — > 0, we can set C2 = 0. This gives the gravitational potential to 
be 

Ci(k)J v/2 (y/wkTj) 5 + 3W 

$w) = ba ; " = 1 ,0 ( 64 ) 

The corresponding density perturbation will be: 

^=-2* g ^72 + (1 + 3u;)V«*»7 (65) 

To understand the nature of the solution, note that du = (a/a)^ 1 oc r\ and kdn — dn/X oc kr\. So the argument of 
the Bessel function is just the ratio (d#/A). From the theory of Bessel functions, we know that for small values of the 
argument J v (x) oc x v is a power law while for large values of the argument it oscillates with a decaying amplitude: 

v j r \ cos{x- (v- 1)tt/4] 

hm J v / 2 {x) K -= — LJ — L ; (66) 

Hence, for modes which are still outside the Hubble radius (k <C i] _1 ), we have a constant amplitude for the potential 
and density contrast: 

$«$i(k); fw-2$ ( (k) (67) 

That is, the perturbation is frozen (except for a decaying mode) at a constant value. On the other hand, for modes 
which are inside the Hubble radius (k 3> r/^ 1 ), the perturbation is rapidly oscillatory (if 10/O). That is the pressure 
is effective at small scales and leads to acoustic oscillations in the medium. 
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A special case of the above is the flat, matter-dominated universe with w = 0. In this case, we need to take the 
w — > limit and the general solution is indeed a constant <!> = ^(k) (plus a decaying mode & decay oc ?7~ 5 which 
diverges as r) — > 0). The corresponding density perturbations is: 

S=-(2 + ^f)^(k) (68) 
which shows that density perturbation is "frozen" at large scales but grows at small scales: 

s = f-2$. i (k) = constant {kq <C 1) 
| -ifc 2 r] 2 $i(k) oc rj 2 oc a (fcry > 1) 

We will use these results later on. 



IV. PERTURBATIONS IN DARK MATTER AND RADIATION 



We shall now move on to the more realistic case of a multi-component universe consisting of radiation and colli- 
sionless dark matter. (For the moment we are ignoring the baryons, which we will study in Sec. IVIfl . It is convenient 
to use y — a/a eq as independent variable rather than the time coordinate. The background expansion of the universe 
described by the function a(t) can be equivalently expressed (in terms of the conformal time rf) as 



PM = a _ x 2, 

PR a cq 

It is also useful to define a critical wave number k c by: 



2x, 



M 



la, 



cq 



1/2 



4(v / 2-l) 2 ??cq 2 -4(V2-l) 2 fc 1 



2 . 

cq > 



l$(Q, m h?)- 1 Mpc 



(70) 
(71) 

k c r\ and y « k c r) 



which essentially sets the comoving scale corresponding to matter-radiation equality. Note that 2x 
in the radiation dominated phase while y = (l/4)(fc c ?7) 2 in the matter dominated phase. 

We now manipulate Eqs. I|52|) . I|55jl. J23), <|57|1 governing the growth of perturbations by essentially eliminating the 
velocity. This leads to the three equations 



2/$' + $ 



lfc 2 



y 



3/c 2 1 



2 1 



1 c 

Om + -Or 
V 



(72) 



(l + yK l + ?^6' m = 3(l + y )<S>ii 



3(2 + 3,) e 
2y kl 



(73) 



4fc 2 



(74) 



for the three unknowns <5>,6 m ,5n. Given suitable initial conditions we can solve these equations to determine the 
growth of perturbations. The initial conditions need to imposed very early on when the modes are much bigger than 
the Hubble radius which corresponds to the y <C 1, k — > limit. In this limit, the equations become: 



1 



Sr 



SL w 3$" 



1 



S' R w 4$" + 2$' 



(75) 



We will take $(yi, k) = $j(fc) as given value, to be determined by the processes that generate the initial perturbations. 
First equation in Eq. Ij75(l shows that we can take 5 R — — 2$^ for yi — > 0. Further Eq. 153JI shows that adiabaticity 
is respected at these scales and we can take 6 m = (3/4)#r = —(3/2)$^;. The exact equation Eq. I|72|l determines $' 
if ($,i5 m ,(5ij) are given. Finally we use the last two equations to set S' m — 3<&',<5' H = 4$', Thus we take the initial 
conditions at some y = yi -C 1 to be: 



$( W ,fe)=#i(fe); 5 R {y h k) = -2®i{k); 8 m {y u k) = -(3/2)$,(fc) 



(76) 



with 5 ' m (y h k) = 3$ ' ( yi , k) ; 5 ' R (y ( , *) = 4$ ' (y h k) . 

Given these initial conditions, it is fairly easy to integrate the equations forward in time and the numerical results 
are shown in Figs [3 El El HI (In the figures k eq is taken to be a eq H eq .) To understand the nature of the evolution, it 
is, however, useful to try out a few analytic approximations to Eqs. (|72|) - l|74|) which is what we will do now. 
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Log 10 (a/a eq ) Log 10 (a/a ecl ) 



FIG. 2: Left PanebThe evolution of gravitational potential $ for 3 different modes. The wavenumber is indicated by the label 
and the epoch at which the mode enters the Hubble radius is indicated by a small arrow. The top most curve is for a mode 
which stays outside the Hubble radius for most of its evolution and is well described by Eq. 1781 . The other two modes show 
the decay of $ after the mode has entered the Hubble radius in the radiation dominated epoch as described by Eq. 11791 . Right 
Panel: Evolution of entropy perturbation (see Eq. 1871 for the definition). The entropy perturbation is essentially zero till the 
mode enters Hubble radius and grows afterwards tracking the dominant energy density perturbation. 



A. Evolution for A 2> dn 



Let us begin by considering very large wavelength modes corresponding to the krj —> limit. In this case adiabaticity 
is respected and we can set Sr ~ (4/3)<5 m . Then Eqs. l(72|) . l(75ll become 



3?/ 



8(1 + y) 



S' R w 4$' 



(77) 



Differentiating the first equation and using the second to eliminate 8 m , we get a second order equation for <E>. Fortu- 
nately, this equation has an exact solution 



1 



16V(l+2/) + 9y 3 + 2y 2 - 8y - 16 



S R « 4$ - 6$, 



(78) 



[There is simple way of determining such an exact solution, which we will describe in Sec. IIVDI ], The initial condition 
on Sr is chosen such that it goes to — 2$j initially. The solution shows that, as long as the mode is bigger than the 
Hubble radius, the potential changes very little; it is constant initially as well as in the final matter dominated phase. 
At late times (y ^S> 1) we see that <!> ~ (9/10)$i so that <& decreases only by a factor (9/10) during the entire evolution 
if k — > is a valid approximation. 



B. Evolution for A <C du in the radiation dominated phase 



When the mode enters Hubble radius in the radiation dominated phase, we can no longer ignore the pressure terms. 
The pressure makes radiation density contrast oscillate and the gravitational potential, driven by this, also oscillates 
with a decay in the overall amplitude. An approximate procedure to describe this phase is to solve the coupled Sr — $ 
system, ignoring 5 m which is sub-dominant and then determine S m using the form of 

When S m is ignored, the problem reduces to the one solved earlier in Eqs <|64[1 . i|fc>5[) with w — 1/3 giving v = 3. 
Since J3/2 can be expressed in terms of trigonometric functions, the solution given by Eq. (|64|) with v = 3, simplifies 
to 

$ = ^ t sin (^) - l v cos (^)] ; P = ^2 ( 79 ) 

Note that as y — > 0, we have $ = $i, $' = 0. This solution shows that once the mode enters the Hubble radius, the 
potential decays in an oscillatory manner. For ly ^> 1, the potential becomes $ w — 3$i(Zy)~ 2 cos(Zy). In the same 



16 



limit, we get from Eq. I|65|l that 

Sr « -|fcV* « « 6$, cos(ly) (80) 

(This is analogous to Eq. (|68|l for the radiation dominated case.) This oscillation is seen clearly in Fig [3] and FigQ] 
(left panel). The amplitude of oscillations is accurately captured by Eq. H80J1 for k = I00k eq mode but not for k — k eq ; 
this is to be expected since the mode is not entering in the radiation dominated phase. 

Let us next consider matter perturbations during this phase. They grow, driven by the gravitational potential 
determined above. When y <C 1, Eq.J73J) becomes: 

C + -C=3$" + ^'-|!$ (81) 

y y K 

The $ is essentially determined by radiation and satisfies Eq. (|61|l ; using this, we can rewrite Eq. (|81|l as 

^(A) = -9(1>' + ^) (82) 

The general solution to the homogeneous part of Eq. (|82|l (obtained by ignoring the right hand side) is (ci + ci In y) ; 
hence the general solution to this equation is 

5 m = ( Cl +c 2 \ny)-9 J ^ J" dyi [$'(yi) + ^ 2 yi*(yi)] (83) 

For j/< 1 the growing mode varies as lnj/ and dominates over the rest; hence we conclude that, matter, driven by $, 
grows logarithmically during the radiation dominated phase for modes which are inside the Hubble radius. 



C. Evolution in the matter dominated phase 



Finally let us consider the matter dominated phase, in which we can ignore the radiation and concentrate on 
Eq. J72J) and Eq. l(75|) . When y > 1 these equations become: 



These have a simple solution which we found earlier (see Eq. (JS5J) 



2k 2 

$ = $oo = const.; S m = -2$^ - ^ $ooy ~ V (85) 



In this limit, the matter perturbations grow linearly with expansion: S m cx y cx a. In fact this is the most dominant 
growth mode in the linear perturbation theory. 




FIG. 3: Evolution of Sr for a mode with k = 100fc cq . The mode remains frozen outside the Hubble radius at (k/k e q) 3 ^ 2 (—Sn) ~ 
(fe/fe e( j) 3/ ' 2 2$ = 2 (in the normalisation used in Fig. and oscillates when it enters the Hubble radius. The oscillations are 
well described by Eq. I8U1 with an amplitude of 6. 
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FIG. 4: Evolution of <5_r for two modes k — fc cq and k = 0.01 fc eq . The modes remain frozen outside the Hubble radius at 
(—Sr) « 2 and oscillates when it enters the Hubble radius. The mode in the right panel stays outside the Hubble radius for 
most part of its evolution and hence changes very little. 



D. An alternative description of matter-radiation system 

Before proceeding further, we will describe an alternative procedure for discussing the perturbations in dark matter 
and radiation, which has some advantages. In the formalism we used above, we used perturbations in the energy 
density of radiation (Sr) and matter (S m ) as the dependent variables. Instead, we now use perturbations in the total 
energy density, S and the perturbations in the entropy per particle, a as the new dependent variables. In terms of 
Sr, S mi these variables are defined as: 

+ yS 

m Pm a /or\ 

a = = = — ; y = — = (86) 

Ptotal PR + Pm l+y PR a cq 

6s\ 36T R 5p m 3 5n R Sn m 

-b R - 6 m = (87) 



s / T R p m 4 n R 




FIG. 5: Evolution of |<5 m | for 3 different modes. The modes are labelled by their wave numbers and the epochs at which they 
enter the Hubble radius are shown by small arrows. All the modes remain frozen when they are outside the Hubble radius and 
grow linearly in the matter dominated phase once they are inside the Hubble radius. The mode that enters the Hubble radius 
in the radiation dominated phase grows logarithmically until y = y cq . These features are well approximated by Eqs. 1831 . 1851 . 
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Given the equations for 5n,5 m , one can obtain the corresponding equations for the new variables (<5, a) by straight 
forward algebra. It is convenient to express them as two coupled equations for $ and a. After some direct but a bit 
tedious algebra, we get: 

V® o/, 2n*/ 3c 2 $ ,fc 2 y ^ 3c 2 a 

y$ + 7777— r + 3(1 + Cg)$ + 777^ r + CsT^TT = 



2(1 + 1/) " 4(1 + y) s fc c 2 l + y 2(1 + y) 



ya' , 2 , , 3c 2 y 2 k 2 _ c 2 y 3 ^ fc ^ 



2 _ (4/3)/* _ 1 A , 3p m ^ _ 1 /, , 3 



2(1 +y) 4(l + y)/c 2 2(l + y)V*c' 

where we have defined 

^=i^ = H 1+ iSJ (90) 

These equations show that the entropy perturbations and gravitational potential (which is directly related to total 
energy density perturbations) act as sources for each other. The coupling between the two arises through the right 
hand sides of Eq. (|88|l and Eq. (|89p. We also see that if we set a = as an initial condition, this is preserved to 
0(k 4 ) and — for long wave length modes — the <I> evolves independent of a. The solutions to the coupled equations 
obtained by numerical integration is shown in Fig. (J2J right panel. The entropy perturbation a ~ till the mode enters 
Hubble radius and grows afterwards tracking either 5r or S m whichever is the dominant energy density perturbation. 
To illustrate the behaviour of <E>, let us consider the adiabatic perturbations at large scales with a w 0, k — ► 0; then 
the gravitational potential satisfies the equation: 

y <p» + + 3(1 + c 2 W + 3c ' $ = 3c > „ (91) 

V + 2(l + y) + ^ i + Cj * + 4(l + y) 2(1 + y) 

which has the two independent solutions: 



both of which diverge as y — ► 0. We need to combine these two solutions to find the general solution, keeping in mind 
that the general solution should be nonsingular and become a constant (say, unity) as y — > 0. This fixes the linear 
combination uniquely: 

■If +-f = — 

HT 1 5 J2 10y 3 



f(y) = —h + 7/2 = T7TT 16^(1 +y) + 9y 3 + 2y 2 -8y- 16 (93) 



Multiplying by $i we get the solution that was found earlier (see Eq. Ij78|0 . Given the form of $ and cr ~ we can 
determine all other quantities. In particular, we get: 

-2(1 + y )d(y<P)/dy + ya 2(1 + y) d 

6r = i+W% "-TTW%^ (y$) (94) 

The corresponding velocity field, which we quote for future reference, is given by: 

2(a/a) ay 

We conclude this section by mentioning another useful result related to Eq. I|88|l . When a as 0, the equation for $ 
can be re-expressed as 

dC 2c 2 k 2 /a 2 p ,„ k 

a± = —J.I ^$«0 for— «1 96 

da i H z p + p ah 



where we have defined: 



C=|^(* + ^+*=^|>° + * (97) 
6 p + p a \ a J P + P k 
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(The i factor arises because of converting a gradient to the k space; of course, when everything is done correctly, all 
physical quantities will be real.) Other equivalent alternative forms for £, which are useful are: 

2 d , *x * H 2 d fa® 



3[1 + w(a)] da y ' a(p + p)dt\H 

For modes which are bigger than the Hubble radius, Eq. (|96|l shows that C is conserved. When £=constant, we can 
integrate Eq. l|9"5)) easily to obtain: 



H 

$ = ci h c 2 



H r da! 

1 



a J H(a') 



(99) 



This is the easiest way to obtain the solution in Eq. I|78|l . 

The conservation law for £ also allows us to understand in a simple manner our previous result that $ only 
deceases by a factor (9/10) when the mode remains bigger than Hubble radius as we evolve the equations from 
y <C 1 to y 3> 1. Let us compare the values of £ early in the radiation dominated phase and late in the matter 
dominated phase. From the first equation in Eq. I|98|l . [using <&' i=s 0] we find that, in the radiation dominated phase, 
£ sa (l/2)$ 4 + $i = (3/2)$j; late in the matter dominated phase, ( sw (2/3)<&/ + <i>/ = (5/3)$/. Hence the conservation 
of C gives $/ = (3/5)(3/2)<i>i = (9/10)$; which was the result obtained earlier. The expression in Eq. I|99|l also works 
at late times in the A dominated or curvature dominated universe. 

One key feature which should be noted in the study of linear perturbation theory is the different amount of growths 
for <E>, Sr and S m . The $ either changes very little or decays; the Sr grows in amplitude only by a factor of few. The 
physical reason, of course, is that the amplitude is frozen at super-Hubble scales and the pressure prevents the growth 
at sub-Hubble scales. In contrast, S m , which is pressureless, grows logarithmically in the radiation dominated era and 
linearly during the matter dominated era. Since the later phase lasts for a factor of 10 in expansion, we get a fair 
amount of growth in S rn . 



V. TRANSFER FUNCTION FOR MATTER PERTURBATIONS 



We now have all the ingredients to evolve the matter perturbation from an initial value 8 = 5i at y = yi -C 1 to the 



current epoch y = ya 



in the matter dominated phase at y 3> 1. Initially, the wavelength of the perturbation 



will be bigger than the Hubble radius and the perturbation will essentially remain frozen. When it enters the Hubble 
radius in the radiation dominated phase, it begins to grow but only logarithmically (see section IIV Bl ) until the 
universe becomes matter dominated. In the final matter dominated phase, the perturbation grows linearly with 
expansion factor. The relation between final and initial perturbation can be obtained by combining these results. 

Usually, one is more interested in the power spectrum Pk(t) and the power per logarithmic band in k— space Ak(t). 
These quantities are defined in terms of Sk(t) through the equations: 



p k {t) = Ai(t) = 



_ k 3 p k (t) 



2tt 2 



(100) 



It is therefore convenient to study the evolution of k 3 ^ 2 Sk since its square will immediately give the power per 
logarithmic band A| in k— space. 

Let us first consider a mode which enters the Hubble radius in the radiation dominated phase at the epoch a en ter- 
From the scaling relation, a en t/k cx t ent oc a 2 nt we find that y en t = (k eq /k). Hence 



k 3 / 2 S m (k,a = l) =— hx(^L) \ k V' 



S ent (k) cx ln( — ) [k 3 / 2 S ent (k)} 

Keq , 



(101) 



MD RD 



at entry 



where two factors — as indicated — gives the growth in radiation (RD) and matter dominated (MD) phases. Let 
us next consider the modes that enter in the matter dominated phase. In this case, a en t/k oc t en f oc a 3 e nt so that 
Vent = {k eq /k) 2 . Hence 



k^ 2 S m (k,a = 1) = — ^ 3 / 2 <Wfc)l « k 2 \e' 2 5 ent {k)] 

Q>ent 



(102) 



MD 



at entry 
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To proceed further, we need to know the k— dependence of the perturbation when it enters the Hubble radius which, 
of course, is related to the mechanism that generates the initial power spectrum. The most natural choice will be 
that all the modes enter the Hubble radius with a constant amplitude at the time of entry. This would imply that the 
physical perturbations are scale invariant at the time of entering the Hubble radius, a possibility that was suggested 
by Zeldovich and Harrison 01 (years before inflation was invented!). We will see later that this is also true for 
perturbations generated by inflation and thus is a reasonable assumption at least in such models. Hence we shall 
assume 



k 3 \S ent (k)\ 2 = k 3 P ent (k) = C = constant, 
Using this we find that the current value of perturbation is given by 



P(k, a = 1) oc S m (k, a = 1) 
The corresponding power per logarithmic band is 

A 2 (k,a = 1) oc k 3 S m (k,a= 1) 



k (for k <C k eq ) 
fc- 3 (lnfc) 2 (for k->k eq ) 



k 4 (for k < k eq ) 
(lnfc) 2 (for k > k eq ) 



(103) 



(104) 



(105) 



The form for P(k) shows that the evolution imprints the scale k eq on the power spectrum even though the initial 
power spectrum is scale invariant. For k < k eq (for large spatial scales), the primordial form of the spectrum is 
preserved and the evolution only increases the amplitude preserving the shape. For k > k eq (for small spatial scales), 
the shape is distorted and in general the power is suppressed in comparison with larger spatial scales. This arises 
because modes with small wavelengths enter the Hubble radius early on and have to wait till the universe becomes 
matter dominated in order to grow in amplitude. This is in contrast to modes with large wavelengths which continue 
to grow. It is this effect which suppresses the power at small wavelengths (for k > k eq ) relative to power at larger 
wavelengths. 



VI. TEMPERATURE ANISOTROPIES OF CMBR 



We shall now apply the formalism we have developed to understand the temperature anisotropics in the cosmic 
microwave background radiation which is probably the most useful application of linear perturbation theory. We 
shall begin by developing the general formulation and the terminology which is used to describe the temperature 
anisotropics. 

Towards every direction in the sky, n = (0, VO we can define a fractional temperature fluctuation A(n) = 
(AT/T)(0,ip). Expanding this quantity in spherical harmonics on the sky plane as well as in terms of the spatial 
Fourier modes, we get the two relations: 

at °° f d 3 k 

A(n) = —(0^)=J2a lm Y lm (B^)= I — -A(k)e^ ni (106) 

where L = t/q — 7y LSS is the distance to the last scattering surface (LSS) from which we are receiving the radiation. The 
last equality allows us to define the expansion coefficients ai m in terms of the temperature fluctuation in the Fourier 
space A(k) . Standard identities of mathematical physics now give 

/H 3 k 
^3 (4tt) i l A(k) 3l {kL) Y lm (k) (107) 

Next, let us consider the angular correlation function of temperature anisotropy, which is given by: 

C(a) = (A(n)A(m)) = £ ^(a ir „4 m ,)^(n)^ m ,(m). (108) 

where the wedges denote an ensemble average. For a Gaussian random field of fluctuations we can express the ensemble 
average as (at m a*, m ,) = CiSu'5 mm '. Using Eq. (|107fl . we get a relation between Ci and A(fc). Given A(fc), the Cj's 
are given by: 

2 f°° 

Ci = - k 2 dk \A(k)\ 2 jftkL) (109) 
7T Jo 



21 



Further, Eq. (|108(l now becomes: 



C(a) = J2 ^r^aPi(cosa) (110) 



Equation 1)110(1 shows that the mean-square value of temperature fluctuations and the quadrupole anisotropy corre- 
sponding to I = 2 are given by 



V T 7 rms 4tt^ V T ) Q An 

These can be explicitly computed if we know A (A;) from the perturbation theory. (The motion of our local group 
through the CMBR leads to a large I — 1 dipole contribution in the temperature anisotropy. In the analysis of CMBR 
anisotropies, this is usually subtracted out. Hence the leading term is the quadrupole with I — 2.) 

It should be noted that, for a given I, the C\ is the average over all m = — 1,0, For a Gaussian 

random field, one can also compute the variance around this mean value. It can be shown that this variance in Ci is 
2Cf /(2l + 1). In other words, there is an intrinsic root-mean-square fluctuation in the observed, mean value of Ci's 
which is of the order of ACi/Ci ss (21 + l) -1 / 2 . It is not possible for any CMBR observations which measures the Cz's 
to reduce its uncertainty below this intrinsic variance — usually called the "cosmic variance". For large values of /, 
the cosmic variance is usually sub-dominant to other observational errors but for low I this is the dominant source of 
uncertainty in the measurement of C;'s. Current WMAP observations are indeed only limited by cosmic variance at 
low- 1. 

As an illustration of the formalism developed above, let us compute the CVs for low I which will be contributed 
essentially by fluctuations at large spatial scales. Since these fluctuations will be dominated by gravitational effects, we 
can ignore the complications arising from baryonic physics and compute these using the formalism we have developed 
earlier. 

We begin by noting that the redshift law of photons in the unperturbed Friedmann universe, v$ — v(a)/a, gets 
modified to the form v$ = v(a)/[a(\ + <E>)] in a perturbed FRW universe. The argument of the Planck spectrum will 
thus scale as 

isq v ( a ) v ( a ) ^( a ) 



T aT (l + $) o(T }[1+($r/4)](1 + *) o(T )[1 + $ + (6 R /4)] 
This is equivalent to a temperature fluctuation of the amount 



(112) 



at large scales. (Note that the observed AT/T is not just (8r/A) as one might have naively imagined.) To proceed 
further, we recall our large scale solution (see Eq. ((75)) ) for the gravitational potential: 



1 



16x7(1 + y) + 9y 3 + 2y 2 - 8y - 16 ; S R - 4$ - 6$, (114) 



I0y 3 

At y = ydec we can take the asymptotic solution $>dec ~ (9/10)$,. Hence we get 



AT 



-5 R + $ 



3 310 1 , N 

= 2$ dec - « 2$ dec - - — $ dec = -$ dec (115) 
A Z y o 



We thus obtain the nice result that the observed temperature fluctuations at very large scales is sim ply related to 
the fluctuations of the gravitational potential at these scales. (For a discussion of the 1/3 factor, see 19] ). Fourier 
transforming this result we get A(k) = (l/3)<i>(k, 77 LSS ) where r] LSS is the conformal time at the last scattering 
surface. (This contribution is called Sachs- Wolfe effect.) It follows from Eq. (|109f) that the contribution to Ci from 
the gravitational potential is 

dk fc 3 |$ fc | 



Q = - I k 2 dk\A(k)\ 2 jUkL) = 1 / **JM_ j ?(feL) (116) 

7T J 7T J k 9 



with 

L = Vo- 7? LSS w Vo w 2(n m H%)- 1 ' 2 » 6000 Q^ 2 h' 1 Mpc (117) 
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For a scale invariant spectrum, fc 3 |$fc| 2 is a constant independent of k. (Earlier on, in Eq. i|103|) we said that scale 
invariant spectrum has fc 3 |5j.| 2 = constant. These statements are equivalent since 6 ~ —2$ at the large scales because 
of Eq. H85J1 with the extra correction term in Eq. 1)85(1 being about 3 x 10~ 4 for k w L^ 1 ,y = ydec-) As we shall see 
later, inflation generates such a perturbation. In this case, it is conventional to introduce a constant amplitude A and 
write: 

fc3|<j) 12 

A| = 1 *' = A 2 = constant (118) 
Substituting this form into Eq. (f 1 1 f>|) and evaluating the integral, we find that 



1(1 + l)Ci _ (A 
2tt 



2 



(119) 



As an application of this result, let us consider the observations of COBE which measured the temperature fluctuations 
for the first time in 1992. This satellite obtained the RMS fluctuations and the quadrupole after smoothing over an 
angular scale of about 9 C ss 10°. Hence the observed values are slightly different from those in Eq. I|lll|) . We have, 
instead, 



Using Eqs. fTT5f> . (fTT^f) we find that 

(£) g «0*4 (%)_"0*LL (121) 

Given these two measurements, one can verify that the fluctuations are consistent with the scale invariant spectrum 
by checking their ratio. Further, the numerical value of the observed (AT/T) can be used to determine the amplitude 
A. One finds that A«3x 1CP 5 which sets the scale of fluctuations in the gravitational potential at the time when 
the perturbation enters the Hubble radius. 

Incidentally, note that the solution Sr = 4$ — 6<J>i corresponds to S m — (3/4)<5r = 3$ — (9/2)$j. At y = ydec, 
taking $ dec = (9/10)$,, we get S m = 3$ dec - (9/2) (10/9) $ dec = -2$ dec . Since (AT/T) obs = (l/3)$ dec we get 
S m = —6(AT/T) b s . This shows that the amplitude of matter perturbations is a factor six larger that the amplitude 
of temperature anisotropy for our adiabatic initial conditions. In several other models, one gets S m = 0(1)(AT/T) o j, s . 
So, to reach a given level of nonlincarity in the matter distribution at later times, these models will require higher 
values of (AT/T) f, s at decoupling. This is one reason for such models to be observationally ruled out. 

There is another useful result which we can obtain from Eq. I|109f) along the same lines as we derived the Sachs- 
Wolfe effect. Whenever fc 3 |A(fc)| 2 is a slowly varying function of k, we can pull out this factor out of the integral and 
evaluate the integral over jf. This will give the result for any slowly varying fc 3 |A(/c)| 2 



l(l + l)Ci fk 3 \A(k)\ 



2?r V 2tt 



(122) 



This is applicable even when different processes contribute to temperature anisotropics as long as they add in quadra- 
ture. While far from accurate, it allows one to estimate the effects rapidly. 



A. CMBR Temperature Anisotropy: More detailed theory 

We shall now work out a more detailed theory of temperature anisotropies of CMBR so that one can understand 
the effects at small scales as well. A convenient starting point is the distribution function for photons with perturbed 
Planckian distribution, which we can write as: 

f( x <*,r 1 ,E,n a ) = ^=f P (^ K y, /p( e )^2[exp( e /T„)-l]- 1 (123) 

The f p (e) is the standard Planck spectrum for energy e and we take e = aE(l + A) -1 to take care of the perturbations. 
In the absence of collisions, the distribution function is conserved along the trajectories of photons so that df/dr) = 0. 
So, in the presence of collisions, we can write the time evolution of the distribution function as 

df _( aE \ ( aE \ 
drj \T+Aj Ip {l + A) 



d\n(aE) dA 
drj drj 



(124) 



r.nll 
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where the right hand side gives the contribution due to collisional terms. Equivalently, in terms of A, the same 
equation takes the form: 



dA d\n(aE) 
drj dr/ 



1 



aE 



Up 



coll 



dA\ 
dr/ J 



(125) 



coll 



To proceed further, we need the expressions for the two terms on the left hand side. First term, on using the standard 
expansion for total derivative, gives: 

dA _ + 9A dx^ + 9A (IB + dA dn^_ ^ Q A + na£) A , UQ s 
dr] di] dx a di] J^i, dr) dr) 

n a zero O(A 2 )=0 

(Note that we are assuming dA/dE — so that the perturbations do not depend on the frequency of the photon.) 
To determine the second term, we note that it vanishes in the unperturbed Friedmann universe and arises essentially 
due to the variation of Both the intrinsic time variation of $ as well as its variation along the photon path will 
contribute, giving: 

dhx(aE) 



dr) 



= -n a d a <S> + <9„$ 



(127) 



(The minus sign arises from the fact that the we have (1 + 2$) in goo but (1 — 2$) in the spatial perturbations. 
Putting all these together, we can bring the evolution equation Eq. I|125fl to the form: 



dA 



— = -n a d a <5> + 0„* + — 



dr) 



dA 



dr) 



(128) 



coll 



Let us next consider the collision term, which can be expressed in the form: 



f dA 
\adrj 



coll 



= -A>.,.A A> 7 ( -S R j A>.,.(vn) 



= N e a T 



-A + -S R + vn 



(129) 



Each of the terms in the right hand side of the first line has a simple interpretation. The first term describes the 
removal of photons from the beam due to Thomson scattering with the electrons while the second term gives the 
scattering contribution into the beam. In a static universe, we expect these two terms to cancel if A = (1/4)5r which 
fixes the relative coefficients of these two terms. The third term is a correction due to the fact that the electrons 
which are scattering the photons are not at rest relative to the cosmic frame. This leads to a Doppler shift which 
is accounted for by the third term. (We denote electron number density by N e rather than n e to avoid notational 
conflict with n a .) 

Formally, Eq. (|128fl is a first order linear differential equation for A. To eliminate the — N e axA term which is linear 
in A in the right hand side, we use the standard integrating factor exp(— r) where 

t( X ) = P di) (aN e a T ) (130) 
Jo 

We can then formally integrate Eq. (|128(l to get: 



A(n) 



We can write 



e- T {-n a d a <$>) 



dxe 
dr) ) 



-r(x) 



-n a <9 a $ + d n $ + aN e a^ 



v-n 



di) 



($e- r ) + (aN e a T <S>)e~ T + (d v $) e - T 



(131) 



(132) 



On integration, the first term gives zero at the lower limit and an unimportant constant (which does not depend on 
n). Using the rest of the terms, we can write Eq. (|131|l in the form: 



A(n) 



dxe T 



2<9„$ + aN e <j„ I $ + \s R + vn 



rno rno / j 

J d X e- T [2d^] + J d X (e- T aN e a T )\§ + -8 R + v>n 



(133) 
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The first term gives the contribution due to the intrinsic time variation of the gravitational potential along the path 
of the photon and is called the integrated Sachs- Wolfe effect. In the second term one can make further simplifications. 
Note that e~ T is essentially unity (optically thin) for z < z rcc and zero (optically thick) for z > z rcc ; on the other 
hand, N e ax is zero for z < z rcc (all the free electrons have disappeared) and is large for z > z roc . Hence the product 
(aN e e~ T ) is sharply peaked at \ — Xrcc (i-e. at z ~ 10 3 with Az ~ 80). Treating this sharply peaked quantity as 
essentially a Dirac delta function (usually called the instantaneous recombination approximation) we can approximate 
the second term in Eq. I|133|) as a contribution occurring just on the LSS: 

A(n) = fc fl + vn + $) +2 f° d X d n ® (134) 

V 4 / LSS Jr, LSS 

In the second term we have put r = oo for r\ < ?7lss an d t — for r\ > ?7lss- 

Once we know Sr, $ and v on the LSS from perturbation theory, we can take a Fourier transform of this result to 
obtain A(fc) and use Eq. I|1U9|) to compute C;. At very large scales the velocity term is sub-dominant and we get back 
the Sachs- Wolfe effect derived earlier in Eq. (|118fl . For understanding the small scale effects, we need to introduce 
baryons into the picture which is our next task. 



B. Description of photon- baryon fluid 

To study the interaction of photons and baryons in the fluid limit, we need to again start from the continuity equation 
and Euler equation. In Fourier space, the continuity equation is same as the one we had before (see Eq. H52[l): 



\)Sr = -k 2 v R + 3$; 



S B = -k 2 v B + 3$ 



(135) 



The Euler equations, however, gets modified; for photons, it becomes: 



vr = + &)-t(v r - v B )\ 



t = N e a T a 



(136) 



The first two terms in the right hand side are exactly the same as the ones in Eq. (|55l) . The last term is analogous to 
a viscous drag force between the photons and baryons which arises because of the non zero relative velocity between 
the two fluids. The coupling is essentially due to Thomson scattering which leads to the factor f . (The notation, and 
the physics, is the same as in Eq. Ijl30|l ). The corresponding Euler equation for the baryons is: 



where 



R: 



a , t(vr — v b ) 

v B = v B + $+-^- '- 

a R 



PB + PB - 3pB ~ 30 n B h 2 ( a 



PR + PR ^PR 



IO- 3 



(137) 



(138) 



Again, the first two terms in the right hand side of Eq. (|137fl are the same as what we had before in Eq. (|55|l . The last 
term has the same interpretation as in the case of Euler equation Eq. (|136|) for photons, except for the factor R. This 
quantity essentially takes care of the inertia of baryons relative to photons. Note that the the conserved momentum 
density of photon-baryon fluid has the form 



(PR +Pr)vr + (pb +Pb)vb ~ (1 + R)(pr+Pr)vr 

which accounts for the extra factor R in Eq. 11371) . 

We can now combine the Eqs. H135(l . Ijl36|) . !|137|) to obtain, to lowest order in (k/f) the equation: 



(139) 
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with 



F 



R 

(1 + R) 



1 

3(1 + J?) 



(141) 
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An exact solution to this equation is difficult to obtain. However, we can try to understand several features by an 
approximate method in which we treat the time variation of R to be small. In that case, we can drop the R terms on 
both sides of the equation. Since we know that the physically relevant temperature fluctuation is A = (1/4)£r + $, 
we can recast the above equation for A as: 

A + fc 2 c 2 A w -k 2 c 2 s R<S> + 2<i> (142) 

Let us further ignore the time variation of all terms (especially $ on the right hand side). Then, the solution is just 
A = — i?$ + Acos(fcc s ?7 LSS ) + -Bsin(/cc s ?7 LSS ). To fix the initial conditions which determine A and B, we recall that 
early on (77 — > 0), we have A — > $/3 (see Eq. I|115[l ) and corresponding velocity should vanish. This gives the solution: 



\Sr + * = y (1 + 3J2) cos(fcc s ?7 LSS ) - ^R; v = (1 + 3R) c s sin(fcc s 7 ?LSS ) (143) 

(One can do a little better by using WKB approximation in which {kc s r]^ ss ) can be replaced by the integral of kc s 
over r\ but it is not very important.) Given this solution, one can proceed as before and compute the C/'s. Adding the 
effects of [$ + (l/4)<5fl] and that of [v • n] in quadrature and noticing that the angular average of ((v • n) 2 ) = (1/3)?; 2 
we can estimate the Ci for scale invariant ( fc 3 |$fc| 2 = 2ir 2 A 2 ) spectrum to be: 



1(1 + l)Ci = 2it 2 A 2 



^ + Q 3 ^ cos(fc*c s ?7 LSS ) - R 



{ l±plc 2 s sin 2 (k*c s ^ ss )\ (144) 



with k* L w / with L = 770 — ^lss — Vo- The key feature is, of course, the maxima and minima which arises from the 
trigonometric functions. The peaks of Ci are determined by the condition A:*c s 77 LSS = lc s r] LSS /r)o = rn:; that is 

W - — f — J = mrV3(l + z dcc f' 2 » 172n (145) 

More precise work gives the first peak at Z pG ak — 200. It is also clear that because of non zero R the peaks are larger 
when the cosine term is negative; that is, the odd peaks corresponding to n — 1, 3, ... have larger amplitudes than the 
even peaks with n = 2,4, .... 

Incidentally, the above approximation is not very good for modes which enter the Hubble radius during the radiation 
dominated phase since <E> does evolve with time (and decays) in the radiation dominated phase. We saw that $ 1=3 
— 3$i(Zy)~ 2 cos(Zy) asymptotically in this phase (see Eq. jSOjl). From Eq. (|50|) we find that during this phase, for 
modes which are inside the Hubble radius, we can take (5r w 6$, cos(iy), so that A « <5r/4 w (3/2)$, cos(iy). On the 
other hand, at very large scale, the amplitude was A = $/3 = (l/3)(9/10)$j = (3/10)$,. Hence the amplitude of the 
modes that enter the horizon during the radiation dominated phase is enhanced by a factor (3/2)(10/9) = 5, relative 
to the large scale amplitude contributed by modes which enter during matter dominated phase. This is essentially 
due to the driving term $ which is nonzero in the radiation dominated phase but zero in the matter dominated phase. 
(In reality, the enhancement is smaller because the relevant modes have k > fc cq rather than k 3> k cq ; see Figs. Eland 

H) 

If this were the whole story, we will see a series of peaks and troughs in the temperature anisotropies as a function of 
angular scale. In reality, however, there are processes which damp out the anisotropies at small angular scales (large 
-I) so that only the first few peaks and troughs are really relevant. We will now discuss two key damping mechanisms 
which are responsible for this. 

The first one is the finite width of the last scattering surface which makes it uncertain from which event we are 
receiving the photons. In general, if V(z) is the probability that the photon was last scattered at redshift z, then we 
can write: 

(%) =[dz{ ^ T ) iHhel f\^V(z). (146) 
\ T J ohs J { scattering was at z J w ' 

From Eq. I|41|l we know that V[z) is a Gaussian with width Az = 80. This corresponds to a length scale 

A; = c (f-) Az- (1 + z dcc ) w ffo' 1 AZ 3/2 ~8 (Slh 2 y 1/2 Mpc. (147) 
\az J qi/2 z / 



over which the temperature fluctuations will be smoothed out. 
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It turns out that there is another effect, which is slightly more important. This arises from the fact that the photon- 
baryon fluid is not tightly coupled and the photons can diffuse through the fluid. This diffusion can be modeled as a 
random walk and the root mean square distance traveled by the photon during this diffusion process will smear the 
temperature anisotropies over that length scale. This photon diffusion length scale can be estimated as follows: 



(Ax) 2 = N 

number o f collisions 



At q 2 At 
= —- = -^1(1) 



q(t) a 2 



(148) 



comovingmean f reepath 

Integrating, we find the mean square distance traveled by the photon to be 



x 2 = 



dt . . _ 3 tdec<l{taec) 



a 2 (t) yw 5 a 2 (t doc ) 

The corresponding proper length scale below which photon diffusion will wipe out temperature anisotropies is: 



(149) 



9diff = a(t d ec)x 



- ^dcc*?(^dcc) 
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35 Mpc 



fl B h 2 
0.02 



(n/4) 
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It turns out that this is the dominant sources of damping of temperature anisotropies at large I rs 10 3 



VII. GENERATION OF INITIAL PERTURBATIONS FROM INFLATION 



In the description of linear perturbation theory given above, we assumed that some small perturbations existed 
in the early universe which are amplified through gravitational instability. To provide a complete picture we need a 
mechanism for generation of these initial perturbations. One such mechanism is provided by inflationary scenario which 
allows for the quantum fluctuations in the field driving the inflation to provide classical energy density perturbations 
at a late epoch. (Originally inflationary scenarios were suggested as pseudo-solutions to certain pseudo-problems; 
that is only of historical interest today and the only reason to take the possibility of an inflationary phase in the early 
universe seriously is because it provides a mechanism for generating these perturbations.) We shall now discuss how 
this can come about. 

The basic assumption in inflationary scenario is that the universe underwent a rapid — nearly exponential - 
expansion for a brief period of time in very early universe. The simplest way of realizing such a phase is to postulate 
the existence of a scalar field with a nearly flat potential. The dynamics of the universe, driven by a scalar field 
source, is described by: 



?(§! =" J "' = 



1 



3A/ 2 ! 



dt 2 dt dd> 



(151) 



where M p i — (8nG) 1 / 2 . If the potential is nearly flat for certain range of <f>, we can introduce the " slow roll-over" 
approximation, under which these equations become: 



H- 



3M 2 ! ' 



3H 



~dt 



(152) 



For this slow roll-over to last for reasonable length of time, we need to assume that the terms ignored in the Eq. I|151|) 
are indeed small. This can be quantified in terms of the parameters: 



<<t>) 



V 



(153) 



which are taken to be small. Typically the inflation ends when this assumption breaks down. If such an inflationary 
phase lasts up to some time t en d then the universe would have undergone an expansion by a factor exp N(t) during 
the interval (t, t en d) where 



a(t) 



Hdt 



1 



A/ 2 , 
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V' 



(154) 
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One usually takes N ~ 65 or so. 

Before proceeding further, we would like to make couple of comments regarding such an inflationary phase. To 
begin with, it is not difficult to obtain exact solutions for a(t) with rapid expansion by tailoring the potential for the 
scalar field. In fact, given any a(t) and thus a H{t) = (a / a) , one can determine a potential V(<f>) for a scalar field 
such that Eq. (|151|l are satisfied (see the first reference in [23 ) • One can verify that, this is done by the choice: 

-2 dH 



I 2 dH~\ r '> ■' r - rl 



(155) 



Given any H(t), these equations give (</>(£), V(t)) and thus implicitly determine the necessary V(<f>). As an example, 
note that a power law inflation, a(t) = aot p (with p 3> 1) is generated by: 

V(4>) = V exp(-^^-) (156) 



P 



pi 



while an exponential of power law 



a{t) oc exp(^), /=4^g. < / < 1, A > (157) 



n<»x( — ) ll-^^H (ir,s, 



can arise from 

V(£\ or I 

,Mpi7 V 6 

Thus generating a rapid expansion in the early universe is trivial if we are willing to postulate scalar fields with tailor 
made potentials. This is often done in the literature. 

The second point to note regarding any inflationary scenarios is that the modes with reasonable size today originated 
from sub-Planck length scales early on. A scale Ao today will be 

X end = \o— = Ao^ » A x ICr 28 (159) 

^0 J end 

at the end of inflation (if inflation took place at GUT scales) and 

Xbegin = KndA^ 1 « A X lO^CA/lO 30 )- 1 (160) 

at the beginning of inflation if the inflation changed the scale factor by A ~ 10 30 . Note that Xbegin < Lp for Ao < 3 
Mpc!! Most structures in the universe today correspond to transplanckian scales at the start of the inflation. It is not 
clear whether we can trust standard physics at early stages of inflation or whether transplanckian effects will lead to 
observable effects j2(J . 

Let us get back to conventional wisdom and consider the evolution of perturbations in a universe which underwent 
exponential inflation. During the inflationary phase the a(t) grows exponentially and hence the wavelength of any 
perturbation will also grow with it. The Hubble radius, on the other hand, will remain constant. It follows that, 
one can have situation in which a given mode has wavelength smaller than the Hubble radius at the beginning of the 
inflation but grows and becomes bigger than the Hubble radius as inflation proceeds. It is conventional to say that 
a perturbation of comoving wavelength Ao "leaves the Hubble radius" when Aoa = dn at some time t = teait(Ao)- 
For t > texu the wavelength of the perturbation is bigger than the Hubble radius. Eventually the inflation ends and 
the universe becomes radiation dominated. Then the wavelength will grow (cx t 1 ' 2 ) slower than the Hubble radius 
(oc t) and will enter the Hubble radius again during t = t enter (Xo). Our first task is to relate the amplitude of the 
perturbation at t = t eX it(Xo) with the perturbation at t = t enter (Ao). 

We know that for modes bigger than Hubble radius, we have the conserved quantity (see Eq. I|97|l 

2 p a ( ■ a \ H ik a . n , 

(=o-J~^ $+-$ +*=_ 5T° + $ (161) 
6 p + p a \ a J P + P K 

At the time of re-entry, the universe is radiation dominated and Centry ~ (2/3)$. On the other hand, during inflation, 
we can write the scalar field as a dominant homogeneous part plus a small, spatially varying fluctuation: cj>(t, x) = 
00 (i) + /(^ x )- Perturbing the equation in Eq. (|151|l for the scalar field, we find that the homogeneous mode <f>Q 
satisfies Eq. Ijl51|l while the perturbation, in Fourier space satisfies: 
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Further, the energy momentum tensor for the scalar field gives [with the "dot" denoting (d/drj) = a(d/dt)\. 
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It is easy to see that $ is negligible at t 
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(164) 



(165) 



(166) 



Thus, given a perturbation of the scalar field //. during inflation, we can compute its value at the time of re-entry, 
which — in turn — can be used to compare with observations. 

Equation (|166fl connects a classical energy density perturbation fk at the time of exit with the corresponding 
quantity <E>fc at the time of re-entry. The next important — and conceptually difficult — question is how we can 
obtain a c-number field fk from a quantum scalar field. There is no simple answer to this question and one possible 
way of doing it is as follows: Let us start with the quantum operator for a scalar field decomposed into the Fourier 
modes with q^t) denoting an infinite set of operators: 



0(t,x) 



ri 3 k 



r«k(t)e* 



(167) 



We choose a quantum state \ip > such the expectation value of quit) vanishes for all non-zero k so that the expec- 
tation value of <p(t, x) gives the homogeneous mode that drives the inflation. The quantum fluctuation around this 
homogeneous part in a quantum state \tp > is given by 

al(t) =< t/>\&{t)\1> > - < > 2 =< M€(t)H > (168) 

It is easy to verify that this fluctuation is just the Fourier transform of the two-point function in this state: 

al(t)= /d 3 x< V#(i,x + y)0(i,y)|V>>e 4k - x . (169) 



Since characterises the quantum fluctuations, it seems reasonable to introduce a c-number field f(t, x) by the 
definition: 



d 3 k 



(170) 



This c-number field will have same c-number power spectrum as the quantum fluctuations. Hence we may take this 
as our definition of an equivalent classical perturbation. (There are more sophisticated ways of getting this result but 
none of them are fundamentally more sound that the elementary definition given above. There is a large literature 
on the subject of quantum to classical transition, especially in the context of gravity; see e.g. [2^]) We now have all 
the ingredients in place. Given the quantum state \ip >, one can explicitly compute Ok and then — using Eq. I|166|) 
with /t = (Tk — obtain the density perturbations at the time of re-entry. 

The next question we need to address is what is \ip >. The free quantum field theory in the Friedmann background 
is identical to the quantum mechanics of a bunch of time dependent harmonic oscillators, each labelled by a wave 
vector k. The action for a free scalar field in the Friedmann background 
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(171) 
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can be thought of as the sum over the actions for an infinite set of harmonic oscillators with mass m = a 6 and 
frequency lo 2 = k 2 /a 2 . (To be precise, one needs to treat the real and imaginary parts of the Fourier transform as 
independent oscillators and restrict the range of k; just pretending that q k is real amounts the same thing.) The 
quantum state of the field is just an infinite product of the quantum state V'klOkj A f° r each of the harmonic oscillators 
and satisfies the Schrodinger equation 

dibir 1 d 2 ib\r 1 9 , 

hr = -2^f + (172) 

If the quantum state V'kfek)^'] of anv given oscillator, labelled by k, is given at some initial time, t' , we can evolve it 
to final time: 

ip [<7k, t}= dq' k K[q k , t; q' k , t'\ij) [q' k , t'} (173) 



where K is known in terms of the solutions to the classical equations of motion and i/j[q' k ,t'} is the initial state. 
There is nothing non-trivial in the mathematics, but the physics is completely unknown. The real problem is that 
unfortunately — in spite of confident assertions in the literature occasionally — we have no clue what ip[g^, if] is. So 
we need to make more assumptions to proceed further. 

One natural choice is the following: It turns out that, Gaussian states of the form 

V>k = A k (t) exp[-B k (t)q£] (174) 

preserve their form under evolution governed by the Schrodinger equation in Eq. I|172f> . Substituting Eq. (|174f) in 
Eq. (|172|l we can determine the ordinary differential equation which governs B k (t). (The A k (t) is trivially fixed by 
normalization.) Simple algebra shows that B k (t) can be expressed in the form 



2 V/k dt 

where fk is the solution to the classical equation of motion 



For the quantum state in Eq. (|174fl . the fluctuations are characterized by 

°l = \(KeB k y 1 = \f k \ 2 (177) 

Since one can take different choices for the solutions of Eq. (|176ll one get different values for a k and different spectra 
for perturbations. Any prediction one makes depends on the choice of mode functions. One possibility is to choose 
the modes so that ip k represents the instantaneous vacuum state of the oscillators at some time t = . (That is Re 
Bh(ti) — (l/2)ui 2 l (ti), say). The final result will then depend on the choice for ti. One can further make an assumption 
that we are interested in the limit of ti — > — oo; that is the quantum state is an instantaneous ground state in the 
infinite past. It is easy to show that this corresponds to choosing the following solution to Eq. I|17tj|) : 

f k = -^=(l + ix)e ix ; x = 4~ (178) 



/2P ' ' Ha 
which is usually called the Bunch-Davies vacuum. For this choice, 

1 A k 2 \ , „ |2| H 2 



\W=M\? + *mp Wl^fcT (179) 

where the second result is at t = t ex it which is what we need to use in Eq. I|ltjtj|) . (Numerical factors of order unity 
cannot be trusted in this computation). We can now determine the amplitude of the perturbation when it re-enters 
the Hubble radius. Eq. (|l(j(jfl gives: 



*k\U ni = [ W r\ W = ¥WS ; k*\* k \i ntty ^l— (iso) 



30 



One sees that the result is scale invariant in the sense that &; 3 |$fc|;? ntry is independent of k. 

It is sometimes claimed in the literature that scale invariant spectrum is a prediction of inflation. This is simply 
wrong. One has to make several other assumptions including an all important choice for the quantum state (about 
which wc know nothing) to obtain scale invariant spectrum. In fact, one can prove that, given any power spectrum 
&{k), one can find a quantum state such that this power spectrum is generated (for an explicit construction, see the 
last reference in |2fT|L So whatever results are obtained by observations can be reconciled with inflationary generation 
of perturbations. 

To conclude the discussion, let us work out the perturbations for one specific case. Let us consider the case of the 
\<j) A model for which 

f/ " •3if^ + y / (<^)=0; y(0)«V o -y0 4 (181) 



Using 



we can write 
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N = 8ttG I r^d<)) = - I -1. - * | M ^— (182) 

/. \—V] 2A > ^ 2 ^ I 0\A,2 



-V(*)=V* (183) 

so that the result in Eq. (|180fl becomes: 

k 3 ' 2 <S> k ~ H 3 A I « X^ 2 N 3 / 2 . (184) 

We do get scale invariant spectrum but the amplitude has a serious problem. If we take N > 50 and note that 
observations require fc 3 / 2 $fc ~ 10 -4 we need to take A < 10~ 15 for getting consistent values. Such a fine tuning of a 
dimensionless coupling constant is fairly ridiculous; but over years inflationists have learnt to successfully forget this 
embarrassment. 

Our formalism can also be used to estimate the deviation of the power spectrum from the scale invariant form. To 
the lowest order we have 

Let us define the deviation from the scale invariant index by (n — 1) = (din A?/dlnfe). Using 

d d d m 2 p V d 

a ~H = 7777 = -7777777 ( 186 ) 



dlnfc da H dcf> 8n V 
one finds that 

1 - n = 6e - 2? 7 (187) 

Thus, as long as e and rj are small we do have n « 1; what is more, given a potential one can estimate e and r\ and 
thus the deviation in — 1). 

Finally, note that the same process can also generate spin-2 perturbations. If we take the normalised gravity wave 
amplitude as h a (, = ^IGirG e a b<^-, the mode function tf> behaves like a scalar field. (The normalisation is dictated by the 
fact that the action for the perturbation should reduce to that of a spin-2 field.) The corresponding power spectrum 
of gravity waves is 



/',„, ( /• ) = = - (—) , n giav (k)h 2 ~ io- 5 (— ) ( i s,si 

2ir z tt \mp J \nip 



Comparing the two results 



a _H^_ ( V 3 \ 2 fV_\ 

scalar ~ {V') 2 ~ \mpV' 2 J ; tonsor ~ \m 2 p J ~ \m% J ( ' 

we get (A tensor /A sca i ar ) 2 w 16ire <C 1. Further, if (1 — n) w 4e (see Eq. (|187f) with e ~ 77) we have the relation 
(A te nsor/A sca i ar ) 2 Ri C(3)(l— n) which connects three quantities, all of which are independently observable in principle. 
If these are actually measured in future it could act as a consistency check of the inflationary paradigm. 



31 



VIII. THE DARK ENERGY 

It is rather frustrating that the only component of the universe which we understand theoretically is the radiation! 
While understanding the baryonic and dark matter components [in particular the values of tts and ^dai] is by no 
means trivial, the issue of dark energy is lot more perplexing, thereby justifying the attention it has received recently. 
In this section we will discuss several aspects of the dark energy problem. 

The key observational feature of dark energy is that — treated as a fluid with a stress tensor T b a = dia (p, — p, —p, —p) 
— it has an equation state p = wp with w < —0.8 at the present epoch. The spatial part g of the geodesic acceleration 
(which measures the relative acceleration of two geodesies in the spacetime) satisfies an exact equation in general 
relativity given by: 

V • g = -4ttG(p + 3p) (190) 

This shows that the source of geodesic acceleration is (p + ip) and not p. As long as (p + 3p) > 0, gravity remains 
attractive while (p + 3p) < can lead to repulsive gravitational effects. In other words, dark energy with sufficiently 
negative pressure will accelerate the expansion of the universe, once it starts dominating over the normal matter. 
This is precisely what is established from the study of high redshift supernova, which can be used to determine the 
expansion rate of the universe in the past [(j . 

The simplest model for a fluid with negative pressure is the cosmological constant (for some recent reviews, see 
\2?fy ) with w = — l,p = —p = constant. If the dark energy is indeed a cosmological constant, then it introduces 
a fundamental length scale in the theory = , related to the constant dark energy density p BE by H\ = 
(8irGp DE /3). In classical general relativity, based on the constants G,c and L\, it is not possible to construct any 
dimensionless combination from these constants. But when one introduces the Planck constant, H, it is possible 
to form the dimensionless combination H^(Gh/c 3 ) = (L 2 P /L\). Observations then require (L 2 P /L\) < 10~ 123 . As 
has been mentioned several times in literature, this will require enormous fine tuning. What is more, in the past, 
the energy density of normal matter and radiation would have been higher while the energy density contributed by 
the cosmological constant does not change. Hence we need to adjust the energy densities of normal matter and 
cosmological constant in the early epoch very carefully so that pa > Pnr around the current epoch. This raises 
the second of the two cosmological constant problems: Why is it that (pa/pnr) = 0(1) at the current phase of the 
universe ? 

Because of these conceptual problems associated with the cosmological constant, people have explored a large variety 
of alternative possibilities. The most popular among them uses a scalar field <fi with a suitably chosen potential V(4>) 
so as to make the vacuum energy vary with time. The hope then is that, one can find a model in which the current 
value can be explained naturally without any fine tuning. A simple form of the source with variable w are scalar fields 
with Lagrangians of different forms, of which we will discuss two possibilities: 

iquin = \da^d a ^ - V(4); itach = -V{<j>)[l - d a cf>d a ^ (191) 

Both these Lagrangians involve one arbitrary function V(4>). The first one, L qu in, which is a natural generalization of 
the Lagrangian for a non-relativistic particle, L = (l/2)q 2 — V(q), is usually called quintessence (for a small sample of 
models, see |24j : there is an extensive and growing literature on scalar field models and more references can be found 
in the reviews in ref.[2j|). When it acts as a source in Friedmann universe, it is characterized by a time dependent 
w(t) with 

P&) = \& + V\ P q (t) = U 2 -y: ™ q = 1 ~ {W,4>2) (192) 
T ' FqK ' T ' q l + (2V/(/) 2 ) 

The structure of the second Lagrangian (which arise in string theory |2^) in Eq. i|191fl can be understood by a 
simple analogy from special relativity. A relativistic particle with (one dimensional) position q(t) and mass m is 
described by the Lagrangian L = — m-y/l — q 2 . It has the energy E = m/ y/l — q 2 and momentum k — mqj \fl — q 2 
which are related by E 2 = k 2 + m 2 . As is well known, this allows the possibility of having massless particles with 
finite energy for which E 2 — k 2 . This is achieved by taking the limit of m — > and q — > 1, while keeping the ratio in 
E — m/yl — q 2 finite. The momentum acquires a life of its own, unconnected with the velocity q, and the energy 
is expressed in terms of the momentum (rather than in terms of q) in the Hamiltonian formulation. We can now 
construct a field theory by upgrading q(t) to a field 4>. Relativistic invariance now requires <f> to depend on both space 
and time [4> = 4>(t, x)] and q 2 to be replaced by di<fid l <fi. It is also possible now to treat the mass parameter m as a 
function of (j>, say, V(<p) thereby obtaining a field theoretic Lagrangian L = —V{4>)y/l — d l (j)di(t>. The Hamiltonian 
structure of this theory is algebraically very similar to the special relativistic example we started with. In particular, 
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the theory allows solutions in which V — > 0, di4>d z (f> — > 1 simultaneously, keeping the energy (density) finite. Such 
solutions will have finite momentum density (analogous to a massless particle with finite momentum k) and energy 
density. Since the solutions can now depend on both space and time (unlike the special relativistic example in which q 
depended only on time), the momentum density can be an arbitrary function of the spatial coordinate. The structure 
of this Lagrangian is similar to those analyzed in a wide class of models called K-essence |26| and provides a rich 
gamut of possibilities in the context of cosmology [27l l2Sj . 

Since the quintessence field (or the tachyonic field) has an undetermined free function V(<j>), it is possible to choose 
this function in order to produce a given H(a). To see this explicitly, let us assume that the universe has two forms 
of energy density with p(a) = pknown(a) + P<j>(a) where Pknown(a) arises from any known forms of source (matter, 
radiation, ...) and prh ja) is due to a scalar field. Let us first consider quintessence. Here, the potential is given 
implicitly by the form |27l |29| . 



6H + 2aH' 



aHQ'l 1 1 1 ' > 1/2 
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where Q(a) = [8irGpk nown (a)/3H 2 (a)] and prime denotes differentiation with respect to a. (The result used in 
Eq. (|155f) is just a special case of this when Q = 0) Given any H(a), Q(a), these equations determine V(a) and </>(a) 
and thus the potential V(<f>). Every quintessence model studied in the literature can be obtained from these equations. 

Similar results exists for the tachyonic scalar field as well |27[. For example, given any H(a), one can construct a 
tachyonic potential V(4>) so that the scalar field is the source for the cosmology. The equations determining V(4>) are 
now given by: 

f da ( aQ' 2aH'\ 1/2 3H 2 , J 2 aH 1 aQ' \ 1/2 , . 

Equations 1)194(1 completely solve the problem. Given any -ff(a), these equations determine V(a) and 4>(a) and thus 
the potential V(cf>). A wide variety of phcnomenological models with time dependent cosmological constant have been 
considered in the literature all of which can be mapped to a scalar field model with a suitable V (</))■ 

While the scalar field models enjoy considerable popularity (one reason being they are easy to construct!) it is very 
doubtful whether they have helped us to understand the nature of the dark energy at any deeper level. These models, 
viewed objectively, suffer from several shortcomings: 

They completely lack predictive power. As explicitly demonstrated above, virtually every form of a(t) can be 
modeled by a suitable "designer" V{4>). 

These models are degenerate in another sense. The previous discussion illustrates that even when w(a) is 
known/specified, it is not possible to proceed further and determine the nature of the scalar field Lagrangian. 
The explicit examples given above show that there are at least two different forms of scalar field Lagrangians 
(corresponding to the quintessence or the tachyonic field) which could lead to the same w(a). (See the second 
paper in rcf. 7] for an explicit example of such a construction.) 

All the scalar field potentials require fine tuning of the parameters in order to be viable. This is obvious in 
the quintessence models in which adding a constant to the potential is the same as invoking a cosmological 
constant. So to make the quintessence models work, we first need to assume the cosmological constant is zero. 
These models, therefore, merely push the cosmological constant problem to another level, making it somebody 
else's problem!. 

By and large, the potentials used in the literature have no natural field theoretical justification. All of them are 
non-renormalisable in the conventional sense and have to be interpreted as a low energy effective potential in 
an ad hoc manner. 

One key difference between cosmological constant and scalar field models is that the latter lead to a w(a) which 
varies with time. If observations have demanded this, or even if observations have ruled out w = —1 at the 
present epoch, then one would have been forced to take alternative models seriously. However, all available 
observations are consistent with cosmological constant (w — —1) and — in fact — the possible variation of w is 
strongly constrained [3(| as shown in Figure El 

While on the topic of observational constraints on w(t), it must be stressed that: (a) There is fair amount of 
tension between WMAP and SN data and one should be very careful about the priors used in these analysis, 
(b) There is no observational evidence for w < — 1. (c) It is likely that more homogeneous, future, data sets 
of SN might show better agreement with WMAP results. (For more details related to these issues, see the last 
reference in |30j . ) 
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FIG. 6: The observational constraints on the variation of dark energy density as a function of redshift from WMAP and SNLS 
data (see [3f|1. The green/hatched region is excluded at 68% confidence limit, red/cross-hatched region at 95% confidence level 
and the blue/solid region at 99% confidence limit. The white region shows the allowed range of variation of dark energy at 
68% confidence limit. 

The observational and theoretical features described above suggests that one should consider cosmological constant 
as the most natural candidate for dark energy. Though it leads to well know fine tuning problems, it also has certain 
attractive features that need to kept in mind. 

• Cosmological constant is the most economical [just one number] and simplest explanation for all the observations. 
We stress that there is absolutely no evidence for variation of dark energy density with redshift, which is 
consistent with the assumption of cosmological constant . 

• Once we invoke the cosmological constant classical gravity will be described by the three constants G, c and 
A = L~^ 2 . It is not possible to obtain a dimensionless quantity from these; so, within classical theory, there 
is no fine tuning issue. Since A(Gft/c 3 ) = (Lp/L\) 2 w 1CP 123 , it is obvious that the cosmological constant is 
telling us something regarding quantum gravity, indicated by the combination Gh. An acid test for any quantum 
gravity model will be its ability to explain this value; needless to say, all the currently available models — strings, 
loops etc. — flunk this test. 

• So, if dark energy is indeed cosmological constant this will be the greatest contribution from cosmology to 
fundamental physics. It will be unfortunate if we miss this chance by invoking some scalar field epicycles! 

In this context, it is worth stressing another peculiar feature of cosmological constant when it is treated as a clue 
to quantum gravity. It is well known that, based on energy scales, the cosmological constant problem is an infra 
red problem par excellence. At the same time, it is a relic of a quantum gravitational effect or principle of unknown 
nature. An analogy |3l| will be helpful to illustrate this point. Suppose one solves the Schrodinger equation for the 
Helium atom for the quantum states of the two electrons tjj(xi,X2)- When the result is compared with observations, 
one will find that only half the states — those in which tp(xi,X2) is antisymmetric under x\ < — > X2 interchange 
- are realized in nature. But the low energy Hamiltonian for electrons in the Helium atom has no information 
about this effect! Here is low energy (IR) effect which is a relic of relativistic quantum field theory (spin-statistics 
theorem) that is totally non perturbative, in the sense that writing corrections to the Helium atom Hamiltonian in 
some (1/c) expansion will not reproduce this result. The current value of cosmological constant could very well be 
related to quantum gravity in a similar way. There must exist a deep principle in quantum gravity which leaves its 
non perturbative trace even in the low energy limit that appears as the cosmological constant . 
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Let us now turn our attention to few of the many attempts to understand the cosmological constant. The choice 
is, of course, dictated by personal bias and is definitely a non-representative sample. A host of other approaches exist 
in literature, some of which can be found in [3^ |. 

A. Gravitational Holography 

One possible way of addressing this issue is to simply eliminate from the gravitational theory those modes which 
couple to cosmological constant. If, for example, we have a theory in which the source of gravity is (p + p) rather 
than (p + 3p) in Eq. (|19t)|l . then cosmological constant will not couple to gravity at all. (The non linear coupling of 
matter with gravity has several subtleties; see eg. |23j.) Unfortunately it is not possible to develop a covariant theory 
of gravity using (p + p) as the source. But we can probably gain some insight from the following considerations. Any 
metric g a b can be expressed in the form g a b = f 2 (x)q a b such that detq = 1 so that det g = f 4 . From the action 
functional for gravity 

it is obvious that the cosmological constant couples only to the conformal factor /. So if we consider a theory of 
gravity in which / 4 = y/— g is kept constant and only q a b is varied, then such a model will be oblivious of direct 
coupling to cosmological constant. If the action (without the A term) is varied, keeping det g = — 1, say, then one 
is lead to a unimodular theory of gravity that has the equations of motion R a b — (l/4:)g a bR = n(T a b — (\/A)g a bT) 
with zero trace on both sides. Using the Bianchi identity, it is now easy to show that this is equivalent to the usual 
theory with an arbitrary cosmological constant. That is, cosmological constant arises as an undetermined integration 
constant in this model 34]. 

The same result arises in another, completely different approach to gravity. In the standard approach to gravity 
one uses the Einstein-Hilbert Lagrangian Leh oc R which has a formal structure Leh ~ R ~ (dg) 2 + d 2 g. If the 
surface term obtained by integrating L sur cx d 2 g is ignored (or, more formally, canceled by an extrinsic curvature 
term) then the Einstein's equations arise from the variation of the bulk term Lbuik oc (dg) 2 which is the non-covariant 
r 2 Lagrangian. There is, however, a remarkable relation between Lbuik and L sur : 

V^gL sur = -d a (ffij^Tp^v^ J (196) 
V d(d agij ) J 

which allows a dual description of gravity using either Lbuik or L sur \ It is possible to obtain the dynamics of gravity 
|35j| from an approach which uses only the surface term of the Hilbert action; we do not need the bulk term at all!. 
This suggests that the true degrees of freedom of gravity for a volume V reside in its boundary dV — a point of 
view that is strongly supported by the study of horizon entropy, which shows that the degrees of freedom hidden 
by a horizon scales as the area and not as the volume. The resulting equations can be cast in a thermodynamic 
form TdS = dE + PdV and the continuum spacetime is like an elastic solid (see e.g. [36^ with Einstein's equations 
providing the macroscopic description. Interestingly, the cosmological constant arises again in this approach as a 
undetermined integration constant but closely related to the 'bulk expansion' of the solid. 

While this is all very interesting, we still need an extra physical principle to fix the value (even the sign) of 
cosmological constant. One possible way of doing this is to interpret the A term in the action as a Lagrange multiplier 
for the proper volume of the spacetime. Then it is reasonable to choose the cosmological constant such that the total 
proper volume of the universe is equal to a specified number. While this will lead to a cosmological constant which 
has the correct order of magnitude, it has several obvious problems. First, the proper four volume of the universe is 
infinite unless we make the spatial sections compact and restrict the range of time integration. Second, this will lead 
to a dark energy density which varies as t~ 2 (corresponding to w = —1/3 ) which is ruled out by observations. 

B. Cosmic Lenz law 

Another possibility which has been attempted in the literature tries to "cancel out" the cosmological constant by 
some process, usually quantum mechanical in origin. One of the simplest ideas will be to ask whether switching on 
a cosmological constant will lead to a vacuum polarization with an effective energy momentum tensor that will tend 
to cancel out the cosmological constant. A less subtle way of doing this is to invoke another scalar field (here we 
go again!) such that it can couple to cosmological constant and reduce its effective value [33 ■ Unfortunately, none 
of this could be made to work properly. By and large, these approaches lead to an energy density which is cither 
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p vv oc L p 4 (where Lp is the Planck length) or to p m oc L A 4 (where La = H A 1 is the Hubble radius associated with 
the cosmological constant ). The first one is too large while the second one is too small! 



C. Geometrical Duality in our Universe 



While the above ideas do not work, it gives us a clue. A universe with two length scales La and Lp will be 
asymptotically De Sitter with a(t) oc exp(i/ La) at late times. There are some curious features in such a universe 
which we will now describe. Given the two length scales Lp and La, one can construct two energy scales p uv = 1/Lp 
and p m = 1/L A in natural units (c = H = 1). There is sufficient amount of justification from different theoretical 
perspectives to treat Lp as the zero point length of spacetime |3Sj . giving a natural interpretation to p vv . The 
second one, p m also has a natural interpretation. The universe which is asymptotically De Sitter has a horizon 
and associated thermodynamics [3^ with a temperature T = and the corresponding thermal energy density 

Pthermai oc T 4 oc 1 / L\ = p m . Thus L p determines the highest possible energy density in the universe while La 
determines the lowest possible energy density in this universe. As the energy density of normal matter drops below 
this value, the thermal ambience of the De Sitter phase will remain constant and provide the irreducible 'vacuum 
noise'. Note that the dark energy density is the the geometric mean p DE = ^/p IR p uv between the two energy densities. 
If we define a dark energy length scale Lp>E such that p DE = l/L% E then Lp,E = \fLpL^ is the geometric mean 
of the two length scales in the universe. (Incidentally, 0.04 mm is macroscopic; it is also pretty close to the 

length scale associated with a neutrino mass of I0~ 2 eV; another intriguing coincidence ?!) 

Using the characteristic length scale of expansion, the Hubble radius dpi = (a/a) -1 , we can distinguish between 
three different phases of such a universe. The first phase is when the universe went through a inflationary expansion 
with dn — constant; the second phase is the radiation/matter dominated phase in which most of the standard 
cosmology operates and dn increases monotonically; the third phase is that of re-inflation (or accelerated expansion) 
governed by the cosmological constant in which dn is again a constant. The first and last phases are time translation 
invariant; that is, t — > t+ constant is an (approximate) invariance for the universe in these two phases. The universe 
satisfies the perfect cosmological principle and is in steady state during these phases! 

In fact, one can easily imagine a scenario in which the two De Sitter phases (first and last) are of arbitrarily long 
duration [4C|. If 17a ~ 0.7, Udm ~ 0.3 the final De Sitter phase does last forever; as regards the inflationary phase, 
nothing prevents it from lasting for arbitrarily long duration. Viewed from this perspective, the in between phase — 
in which most of the 'interesting' cosmological phenomena occur — is of negligible measure in the span of time. It 
merely connects two steady state phases of the universe. 

While the two De Sitter phases can last forever in principle, there is a natural cut off length scale in both of them 
which makes the region of physical relevance to be finite [4(J]. Let us first discuss the case of re-inflation in the late 
universe. As the universe grows exponentially in the phase 3, the wavelength of CMBR photons are being redshifted 
rapidly. When the temperature of the CMBR radiation drops below the De Sitter temperature (which happens when 
the wavelength of the typical CMBR photon is stretched to the La-) the universe will be essentially dominated by 
the vacuum thermal noise |39j due to the horizon in the De Sitter phase. This happens when the expansion factor is 
a = ap determined by the equation To(ao/ap) = (l/2irL\). Let a = a\ be the epoch at which cosmological constant 
started dominating over matter, so that (a\/ao) 3 — (£Idai /^a)- Then we find that the dynamic range of the phase 
3 is 



Interestingly enough, one can also impose a similar bound on the physically relevant duration of inflation. We know 
that the quantum fluctuations, generated during this inflationary phase, could act as seeds of structure formation 
in the universe. Consider a perturbation at some given wavelength scale which is stretched with the expansion of 
the universe as A oc a(t). During the inflationary phase, the Hubble radius remains constant while the wavelength 
increases, so that the perturbation will 'exit' the Hubble radius at some time. In the radiation dominated phase, 
the Hubble radius dn oc t oc a 2 grows faster than the wavelength A oc a(t). Hence, normally, the perturbation will 
're-enter' the Hubble radius at some time. If there was no re-inflation, all wavelengths will re-enter the Hubble radius 
sooner or later. But if the universe undergoes re-inflation, then the Hubble radius 'flattens out' at late times and 
some of the perturbations will never reenter the Hubble radius ! If we use the criterion that we need the perturbation 
to reenter the Hubble radius, we get a natural bound on the duration of inflation which is of direct astrophysical 
relevance. Consider a perturbation which leaves the Hubble radius (H^ 1 ) during the inflationary epoch at a = cti- It 
will grow to the size H~ 1 L (a/ai) at a later epoch. We want to determine ai such that this length scale grows to L\ 
just when the dark energy starts dominating over matter; that is at the epoch a = a\ = ao(Qp,M /^a) 1 ^ 3 • This gives 
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H~^(a\/ai) — £a so that a, = {H~^ j 'L\){Q.dm /^A) 1 ^ 3 ao- On the other hand, the inflation ends at a = a enc i where 
a end/o-o = ^o/^reheat where T^eat is the temperature to which the universe has been reheated at the end of inflation. 
Using these two results we can determine the dynamic range of this phase 1 to be 

a -T = (w 3 ^) (w^) = 2 r OF/0A 2-i * 1025 (198) 

a i V J reheat / \ "-DM/ 2tT1 reheat-"i n 

where we have used the fact that, for a GUTs scale inflation with Eqjjt — 10 14 GeV, T re heat = EcuTjPin = Eq UT 
we have 27ri? i I T re heat = (37T /2) 1 / 2 {Ep / Equt) ~ 10 5 . If we consider a quantum gravitational, Planck scale, inflation 
with 27T-ff "^Tj-oheat = the ranges in Eq. I|197|) and Eq. (|198fl are approximately equal. 

This fact is definitely telling us something regarding the duality between Planck scale and Hubble scale or between 
the infrared and ultraviolet limits of the theory. The mystery is compounded by the fact the asymptotic De Sitter 
phase has an observer dependent horizon and related thermal properties |3Sl | . Recently, it has been shown — in a scries 
of papers, see ref . |35j — that it is possible to obtain classical relativity from purely thermodynamic considerations. 
It is difficult to imagine that these features are unconnected and accidental; at the same time, it is difficult to prove 
a definite connection between these ideas and the cosmological constant. 



D. Gravity as detector of the vacuum energy 



Finally, we will describe an idea which does lead to the correct value of cosmological constant. The conventional 
discussion of the relation between cosmological constant and vacuum energy density is based on evaluating the zero 
point energy of quantum fields with an ultraviolet cutoff and using the result as a source of gravity. Any reasonable 
cutoff will lead to a vacuum energy density p V ac which is unacceptably high. This argument, however, is too simplistic 
since the zero point energy — obtained by summing over the (l/2)hujk — has no observable consequence in any 
other phenomena and can be subtracted out by redefining the Hamiltonian. The observed non trivial features of the 
vacuum state of QED, for example, arise from the fluctuations (or modifications) of this vacuum energy rather than 
the vacuum energy itself. This was, in fact, known fairly early in the history of cosmological constant problem and 
is stressed by Zeldovich |4l| who explicitly calculated one possible contribution to fluctuations after subtracting away 
the mean value. This suggests that we should consider the fluctuations in the vacuum energy density in addressing 
the cosmological constant problem. 

If the vacuum probed by the gravity can readjust to take away the bulk energy density p uv ~ Lp 4 , quantum 
fluctuations can generate the observed value pde- One of the simplest models |43 which achieves this uses the 
fact that, in the semi-classical limit, the wave function describing the universe of proper four-volume V will vary 
as ^ (x exp(— iAq) oc exp[— i(A e gV / Lp)]. If we treat (A/Lp, V) as conjugate variables then uncertainty principle 
suggests AA « Lp/AV. If the four volume is built out of Planck scale substructures, giving V = NLp, then the 
Poisson fluctuations will lead to AV « y/VL p giving AA = Lp/AV ss 1/y/V ~ Hq. (This idea can be a more 
quantitative; see |42^)- 

Similar viewpoint arises, more rigorously, when we study the question of detecting the energy density using gravi- 
tational field as a probe. Recall that an Unruh-DeWitt detector with a local coupling Lj = M{t)4>[x{t)\ to the field 
cf> actually responds to (Q\<j)(x)4>(y)\Q) rather than to the field itself |4^. Similarly, one can use the gravitational field 
as a natural "detector" of energy momentum tensor T a b with the standard coupling L — nh a i,T ab . Such a model 
was analysed in detail in ref. [44J and it was shown that the gravitational field responds to the two point function 
(0\T a b(%)T c d(y)\0) . In fact, it is essentially this fluctuations in the energy density which is computed in the inflationary 
models (see Eq. (117011 s ) as the seed source for gravitational field, as stressed in ref. All these suggest treating the 
energy fluctuations as the physical quantity "detected" by gravity, when one needs to incorporate quantum effects. If 
the cosmological constant arises due to the energy density of the vacuum, then one needs to understand the structure 
of the quantum vacuum at cosmological scales. Quantum theory, especially the paradigm of renormalization group 
has taught us that the energy density — and even the concept of the vacuum state — depends on the scale at which 
it is probed. The vacuum state which we use to study the lattice vibrations in a solid, say, is not the same as vacuum 
state of the QED. 

In fact, it seems inevitable that in a universe with two length scale La, Lp, the vacuum fluctuations will contribute 
an energy density of the correct order of magnitude p DE = ^/p IR p uv . The hierarchy of energy scales in such a universe, 
as detected by the gravitational field has 40, 45] the pattern 



Pvac ~ L% + L p (lI) + L p {l a ) + 



(199) 
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The first term is the bulk energy density which needs to be renormalized away (by a process which we do not 
understand at present); the third term is just the thermal energy density of the De Sitter vacuum state; what is 
interesting is that quantum fluctuations in the matter fields inevitably generate the second term. 

The key new ingredient arises from the fact that the properties of the vacuum state depends on the scale at which it 
is probed and it is not appropriate to ask questions without specifying this scale. If the spacetime has a cosmological 
horizon which blocks information, the natural scale is provided by the size of the horizon, La, and we should use 
observables defined within the accessible region. The operator H{< La), corresponding to the total energy inside a 
region bounded by a cosmological horizon, will exhibit fluctuations AE since vacuum state is not an eigenstate of 
this operator. The corresponding fluctuations in the energy density, Ap oc (AE)/L A — f(Lp,L\) will now depend 
on both the ultraviolet cutoff Lp as well as L\. To obtain Ap vac oc AE/L\ which scales as (LpLa)~ 2 we need to 
have (AE) 2 oc Lp 4 L^; that is, the square of the energy fluctuations should scale as the surfa ce a rea of the bounding 
surface which is provided by the cosmic horizon. Remarkably enough, a rigorous calculation |4o| of the dispersion in 
the energy shows that for La ^> Lp, the final result indeed has the scaling 

(A£) 2 = Cl ^| (200) 

where the constant c\ depends on the manner in which ultra violet cutoff is imposed. Similar calculations have been 
done (with a completely different motivation, in the context of entanglement entropy) by several people and it is 
known that the area scaling found in Eq. I|200l) . proportional to L\, is a generic feature |4y|- F° r a simple exponential 
UV-cutoff, ci = (1/307T 2 ) but cannot be computed reliably without knowing the full theory. We thus find that the 
fluctuations in the energy density of the vacuum in a sphere of radius L\ is given by 

AE H 2 
Apvac = -js- cx L p 2 Ll 2 ex (201) 

The numerical coefficient will depend on c\ as well as the precise nature of infrared cutoff radius (like whether it is La 
or La/27t etc.). It would be pretentious to cook up the factors to obtain the observed value for dark energy density. 
But it is a fact of life that a fluctuation of magnitude Ap vac ~ H\/G will exist in the energy density inside a sphere 
of radius 1 if Planck length is the UV cut off. One cannot get away from it. On the other hand, observations 
suggest that there is a p vac of similar magnitude in the universe. It seems natural to identify the two, after subtracting 
out the mean value by hand. Our approach explains why there is a surviving cosmological constant which satisfies 
P-de = \/PmPvv which — in our opinion — is the problem. 
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